Controlling and measuring quantum transport of heat in trapped-ion crystals 
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Measuring heat flow through nanoscale systems poses formidable practical difficulties as there is no 'ampere 
meter' for heat. We propose to overcome this problem by realizing heat transport through a chain of trapped 
ions. Laser cooling the chain edges to different temperatures induces a current of local vibrations (vibrons). We 
show how to efficiently control and measure this current, including fluctuations, by coupling vibrons to internal 
ion states. This demonstrates that ion crystals provide a suitable platform for studying quantum transport, 
e.g., through thermal analogues of quantum wires and quantum dots. Notably, ion crystals may give access to 
measurements of the elusive large fluctuations of bosonic currents and the onset of Fourier's law. These results 
are supported by numerical simulations for a realistic implementation with specific ions and system parameters. 



In view of the rapid development of nanoscale technolo- 
gies [1], understanding charge and heat transport at the micro- 
scopic level has become a central topic of current research. 
As already shown for fermions [2], charge transport at the 
nanoscale is typically governed by quantum effects. For 
bosons, the transport of heat by phonons has been predicted 
to have analogous properties [3]. However, experiments are 
considerably more challenging in this case, as there is no de- 
vice capable of measuring local heat currents [3]. Moreover, 
the heat reservoirs and temperature probes required to study 
heat transport, usually entail spurious interface effects. Within 
these restrictions, most experimental efforts have focused on 
detecting temperature profiles [4] in different devices [1, 5, 6]. 

In this Letter, we show that trapped-ion crystals are a 
promising platform to overcome these limitations. By exploit- 
ing laser-induced couplings between transverse quantized vi- 
brations (vibrons) and internal degrees of freedom (spins) of 
the ions, we show how to control and measure the heat current 
across an ion chain [Fig. 1(a)]. We introduce a quantum trans- 
port toolbox containing all functionalities required for treating 
heat currents on the same footing as electrical currents. To 
be more precise, ions at the edges of the crystal are Doppler- 
cooled to different temperatures, and act as reservoirs that sus- 
tain a non-equilibrium heat flow through the bulk of the crys- 
tal [Fig. 1(b)]. To probe this transport, we propose to map the 
information of the vibron number or heat current (including 
their fluctuations) onto the spins, which can be measured via 
spin-dependent fluorescence [7]. For conventional materials, 
the quest for such a probing device is formidable, as it would 
require measuring the heat flow through a single atom. 

We demonstrate the versatility of our toolbox with two ex- 
amples: a thermal quantum wire (TQW), and a thermal quan- 
tum dot (TQD). In the TQW, we study the onset of Fourier's 
law [8] in a mesoscopic system, which predicts a temperature 
gradient across a heat conductor. This law requires the tran- 
sition from ballistic to diffusive transport, which can be in- 
duced by either /) dephasing through noisy modulations of the 
trap frequencies [9], or //) disorder due to an engineered spin- 
vibron coupling [10]. The TQD highlights the differences of 
bosonic and fermionic transport [11], captured by the statistics 
of the fluctuations [12]. Building on the idea of laser-assisted 
tunneling [13, 14], we show how to measure cuiTent fluctua- 
tions. Moreover, the TQD can be operated as a switch for the 
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Figure 1. Heat transport toolbox: (a) (upper panel) Mixed-species 
(blue= T, red= a, green = K) ion crystal in a linear Paul trap. Our 
ideas also apply to surface trap arrays, (lower panel) The spin and 
vibrational degrees of freedom of the ions are indicated with ar- 
rows and wells, respectively. We show the laser arrangements L|™ 
(IJ^ ,1}^) to control the incoherent (coherent) vibrational dynamics 
of the ions, where Jg,p represents the vibron tunneling, (b) (upper 
panel) Scheme for a thermal quantum wire (TQW) connected to two 
reservoirs held at different temperatures, (lower panel) A sufficiently 
strong laser cooling with strengths A^ > A^ , allows to treat the edge 
T-ions as heat reservoirs, whereas the bulk CT, K" ions act as the TQW. 



heat current, a first step towards a single-spin heat transistor. 

We note that the control of heat in trapped ions is a topic of 
increasing interest. InRef. [15], the propagation of vibrational 
excitations along an ion chain has been assessed theoretically. 
The use of single trapped ions as heat engines has been pro- 
posed in [16]. Finally, and more relevant to the present topic, 
the thermalization of a sympathetically-cooled ion chain has 
been studied in [17] by Langevin dynamics [18]. We believe 
that the toolbox hereby introduced, which is based on a thor- 
ough first-principle derivation, will be genuinely useful for the 
development of heat transport experiments with trapped ions. 

Model- We consider a linear Coulomb crystal with three 
types of ions [Fig. 1(a)]. The quanta of the collective crystal 
vibrations (i.e. the phonons) are customary used as a data bus 
for quantum information processing [19]. In this work, how- 
ever, we are interested in the individual transverse oscillations 



of the ions, whose quanta will be called vibrons, and are re- 
sponsible for a local electric dipole. As demonstrated in recent 
experiments [20], the interaction between these dipoles leads 
to a tight-binding model for the transverse vibrations {h = \) 
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where the bosonic operators a- [a- ) create (annihilate) the lo- 
cal vibrons, and we have used latin indices to label the lattice 
sites /, J e {1 • • -A^}, and greek sub-indices to label the species 
a,/3 G {(7,T, K"} [Fig. 1(a)]. According to Eq. (1), the ion 
crystal is a natural playground for bosonic lattice models [21], 
where the bosonic particles correspond to the vibrons, the lat- 
tice is determined by the underlying crystal structure, and the 
trapping and dipole-dipole couplings are responsible for the 
on-site energies a),^, and the long-range tunnelings Jjajn [22]. 

In addition to the vibrations, we exploit two selected atomic 
levels of each ion, hereafter referred to as spins \si^) E 
{|t/„),i;,a>}- The spin Hamiltonian is //f = ^I,-„<(J4, 
where af^ = |t/a)(t(al - \iia){iial and the atomic ti-ansition 
is characterized by the frequency co" and the linewidth Fa- 

Let us now augment the vibron and spin dynamics given 
by ii/tb + Ha with additional incoherent and coherent contribu- 
tions to obtain the tools for studying quantum transport: 

(i) For the incoherent dynamics, we employ a laser L™ 
forming a standing wave along the direction of the vibrons. 
This drives a dipole-allowed transition of the a-spins, simulta- 
neously modifying the number of vibrational excitations. For 
a fast decay of the a-spins [22], the interplay of these two 
processes yields an effective dissipation of the vibrons 

^;«(M) -i^[A+,oL,o,J(M) + i^[Ar„,«,v<KM), (2) 

where ^[A,C>i,C>2](») = A(C>i • O2 - 020i») +H.c. is a 
generic super-operator throughout this work. The local heat- 
ing (cooling) strengths A^ (A^ ), which depend on the spec- 
tral functions of the couplings [23], can be tuned by control- 
ling the laser parameters in the Doppler regime [22, 26]. 

(//) For the coherent dynamics, we apply a spin-dependent 
traveling wave L'^, which consists of a pair of non- 
copropagating laser beams. The wave originates from two- 
photon processes, whereby the spin state is virtually ex- 
cited by absorbing a photon from a laser beam, and sub- 
sequently de-excited by photon emission into the remaining 
laser beam [22]. This results in a spin-vibron interaction 

Hi'i{t) = \{^(0++^(0^al)co^{Vat-(pa)a]^a.^, (3) 

where the parameters Ao);^, Va, ^a are fully controllable [29]. 
The Hamiltonian and Lindblad terms in Eqs. (1) to (3) form 
our Liouvillian toolbox for studying quantum transport in ion 
crystals: the driven dissipative spin-vibron model 
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Here, the sums over € and 3 comprise ions subjected to either 
additional coherent or incoherent effects, respectively. Let us 



remark that our toolbox avoids using single-ion laser address- 
ing by exploiting different species a for each functionality. 

A key function of our toolbox is the possibility of designing 
reservoirs. Ideally, a reservoir in quantum transport theory is 
a system capable of supplying and absorbing carriers without 
changing its state. In our case, this can be achieved by using a 
red-detuned laser L™ such that the dissipation in Eq. (2) im- 
plements a cooling process at the rate Yj^ = Re{(A|^ )* — A^ }. 
For sufficiently strong laser intensities, the cooling dominates 
over the tunneling 7;^ ^ JiaJa [22]. Thus, the ion does not 
have enough time to exchange vibrons with the rest of the 
chain before thermalizing. More importantly, the constantly 
cooled ions remain in a thermal state u* with average vibron 

number (equivalently temperature), n,^ — Re{A^ }/Yia, and 
thus provide an accurate realization of a vibronic reservoir 

Thermal quantum wire- We now consider how to design a 
thermal quantum wire (TQW) with our toolbox. We choose 
certain ion species [22], such that dissipation only occurs for 
the T-ions Ft > FctiFic (i.e. CT, ff e £ and T G 3). By plac- 
ing the T-ions at the edges of the chain, the dissipator (2) 
realizes the starting point of different transport studies [31]. 
Considering the strong-cooling regime with h\^> Hn^, the T- 
ions resemble a vibronic battery. The left reservoir constantly 
supplies vibrons in the attempt to equilibrate with the TQW 
by increasing its vibrational excitations. Conversely, the right 
reservoir keeps absorbing vibrons in its effort to lower the vi- 
brational excitations in the TQW. When combined, the reser- 
voirs sustain a flow of heat along the TQW [Fig. 1(b)]. 

To start with, we assess how the TQW thermalizes due to 
the contact with the reservoks on a time-scale much longer 
than yf' ' . Since the vibron tunnelings are much smaller than 
the edge cooling rates, the edge vibrons are "integrated out" 
by projecting onto the state jx -^ ji^^ (g) /Xbuik ^ jU* . We obtain 
a dissipative spin-vibron model for the reduced density matrix 
of the bulk of the wire, namely (9r/ibuik = -S^ddsv (Mbuik), where 
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Here, i/nb is identical to the tight-binding model (1) with 
renormalized tunnelings and on-site energies [22]. The more 
relevant part for the thermalization of the bulk ions is the dis- 
sipation due to tunneling exchange with the reservoirs. This 
is similar to Eq. (2) but extended to all possible pairs of ions 
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where we introduced non-local cooling and heating strengths 
K ; = Y.I M' ; («£. + 1), and At . = \n irf^ . «. , 
respectively. These parameters depend on the sequential 
bulk-edge-bulk tunneling processes via the coupling matrix 
^t-jp ^ ^^Jia.itPiA(^ia)JerJpi whcrc we have introduced a 
Lorentzian density of states for the laser-cooled T-ions 
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The process of thermalization of the TQW becomes clear 
from this derivation: Tunneling of vibrons from the bulk of 



the wire into the reservoirs, and back, introduces an effective 
dissipation responsible for its thermalization. The cubic decay 
of the dipole-dipole tunnehngs with distance suggests that the 
thermaHzation mainly occurs by exchange of vibrons between 
the reservoirs and the neighboring bulk ions. Moreover, since 
Ji^ » Jig_j„ the broad density of states Pf\.(e) does not allow 
to resolve the different vibrational levels of the TQW. We can 
thus predict a homogeneous steady-state vibron occupation 
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where we have defined the system-reservoir couplings as Fl = 



r,J, ,\ ,„ ,> , and the reservoir mean vi- 
: «i^ , «R = nftf^. Similar arguments apply to 



Yr,^ r, and Fr = 

bron numbers hl 

the vibron current, defined through (9,«,„ 

turns out to be independent of the length of the TQW 
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Numerical solutions of the full dissipative dynamics in 
Eq. (5) confirms our predictions [22]. These results violate 
Fourier's law of thermal conduction [8], which predicts: (i) A 
linear temperature gradient, which is equivalent to (n,„)FL = 
"L + (mr — nh)ia/N. (ii) A heat current that is inversely pro- 
portional to the length of the wire (/"'^)fl °^ («l ^ nRj/N. 
The violation of these predictions is not a surprise, since 
Fourier's law only applies to diffusive processes. In our case, 
however, Eqs. (6)-(7) describe the ballistic transport of vi- 
brons through the TQW, and are completely analogous to the 
equivalent expressions for electronic transport [32]. 

We consider two phase-breaking processes for a ballistic- 
diffusive crossover: dephasing [33] and disorder [34]. One 
should not expect a perfect recovery of Fourier's law for the 
small crystals realized in trapped-ion laboratories. However, 
we show that diffusion leads to clear signatures in the vibron 
occupation measurable in state-of-the-art experiments. 

To engineer dephasing, a possibility is to modulate the trap 
frequencies with a noisy voltage [9]. We model such noise as 
a dynamic fluctuation of the on-site energies in Eq. (1), a),-^ — > 
a),„ + 5a),„(f), where 5cOii^{t) is a random process [22]. In a 
Born-Markov approximation, this leads to an additional term 
in Eq. (4), .i^ddsv ~^ -^ddsv + ^d> where we have introduced 
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such that Fd is the dephasing rate, and ^c accounts for a noise 
correlation length [35]. In Fig. 2(a), we show the homoge- 
neous distribution of vibrons along the TQW in the absence 
of dephasing. When dephasing with ^i^<^ Lis switched on, 
we observe the onset of a linear gradient along the microtrap 
array [upper panel of Fig. 2(b)], which pinpoints the presence 
of diffusive transport. For the linear Paul trap, the inhomo- 
geneity of the crystal modifies the linear gradient leading to 
an anomalous Fourier's law [lower panel of Fig. 2(b)]. 

Let us now turn to the problem of disorder. The scatter- 
ing of vibrons by impurities can be modeled by modifying the 



on-site energies of Eq. (1), (B,-^ -^ C0,„ + A(Oa, where AOa is 
a static random variable. To obtain such disorder, we apply 
a strong static spin-vibron coupling (3) with parameters Va = 
0, (pa = 0, Aco^ = 0, and Aco^ 7^ 0, such that the vibrons cross- 
ing the TQW experience an inhomogeneous landscape of on- 
site energies depending on the spin configuration [10]. When 
each bulk spin is initialized in |+,a) — (|tia) + lifa))/v^' '^^e 
tight-binding model (1) becomes stochastic [22], i/tb —^ H^tb 
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Here, the on-site energies for the ions a.ja G 2! are binary 
random variables sampling e,^^ e {(Oa — jAco^ ,(x>a + j^(Oa} 
with a flat probability distribution p(e,„) = 5 inherited from 
the quantum parallelism. This randomness leads to Anderson 
localization, whereby normal modes display a finite localiza- 
tion length ^loc [36]. For the small ion crystals of length L 
here considered, the modes with ^loc ^ L contribute ballisti- 
cally, those with £,ic,c ^ L introduce diffusion, and £,]oc ^ L do 
not contribute to the transport. Hence, one can expect that the 
heat transport will be much richer in the disordered case. In 
the upper (lower) panel of Fig. 2(c), we display the disorder- 
averaged distribution of vibrons along the TQW of a microtrap 
array (linear Paul trap). Far away from the edges, we recog- 
nize again the linear gradient representative of Fourier's law. 
To distinguish between ballistic and diffusive transport we 
require a measuring tool for the vibron occupations. We sug- 
gest a general scheme for mapping the mean value of any vi- 
bron operator ((9,„)ss, and its fluctuation spectrum 



5o,„oJft))=|df(O,„(f)O,,(0)),,e- 



Oia^O„-{di, 



onto the spin coherences, while disturbing the vibron states 
minimally. This idea goes along the lines of exploiting a 
qubit to measure phase fluctuations in a Bose-Einstein con- 
densate [37] or full counting statistics of charge transport [38]. 
By engineering a spin-vibron coupling of the form H^ — 
^U h^oOi^of , where Xq is a weak coupling constant, the 
mean value and the fluctuations of the vibronic operator Ot^ 
can be measured by Ramsey-type interferometry [22]. If the 
probe consists of a single K"-ion [39] initialized by a 7r/2-pulse 
in the spin state |+,,^), the ic-spin acquires phase information 
about the steady-state vibron observable. In order to recover 
this information, we perform another ;r/2-pulse, and measure 
the probability of observing the K"-ion in the spin-down state. 
This is equivalent to the measurement of the spin coherences 



{af^t)) = cos(Ao(0,J.sOe"^^'''^'''-''-'°'^'. 



(8) 



Therefore, the period (decay) of these spin oscillations yields 
the mean value (zero-frequency fluctuations) of the vibron op- 
erator {Oi^)ss (Soj O;-(0))- Due to the excellent accuracies 
of projective spin measurements with ion traps, which ex- 
ceed 99% with millisecond detection times [19], probing the 
steady-state vibrons with this method promises to be very effi- 
cient. Specifically, the scheme allows us to measure the mean 
vibron number Oi^ = rtj^. The required Ramsey -type coupling 
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Figure 2. Fourier's law: Steady-state of a chain with N = 20 ions (upper panels: microtrap array, lower panels: linear Paul trap). We display 
the vibron distribution in (a) the ballistic regime, which is in agreement with the numerical results of [17], (b) the dephasing-induced diffusive 
regime, and (c) the disorder-induced diffusive regime. 



is obtained directly by switching on a weak static spin-vibron 
coupling (3) with parameters v^ — Q,(Pk — 0,Aa)+ = 0, and 
Aco^ 7^ 0. Then, the dynamics of the probe spin gives di- 
rectly the mean vibron number (n;^)sss and its fluctuations 
Srii n, (0). Let us note that, if the probe consists of several 
Jf-ions, we get information about density-density correlators. 

Thermal quantum dot and single-spin heat switch.- For re- 
alizing a thermal quantum dot (TQD), we consider a single 
K"-ion at position Pk: in the bulk of the crystal. We use the 
remaining a-ions as thermal contacts by employing a strong 
static spin-phonon coupling (3) with parameters Va =0,(pa = 
OjAOfjr = 0, and Aco^ ^ 0. When preparing the spins in 
\Wo) = lltj • • • ia)\<l>K) Ita • • • tcj), there is a large shift of the 
on-site energies at p^, inhibiting tunneling across the junction. 
Hence, the two halves of the chain thermalize independently 
(n,-„)ss = "L for ia <Pk, and (n;^)ss = "r, for ic > Pk, leading 
to an analogue of a quantum dot between two leads. The Li- 
ouvillian is .^^^^^ ^^l + -^lkr + ^r, where JfL/R describe 
the uncoupled halves (5), and ^licr describes the TQD. 

So far, the transport through the TQD is blocked by the 
large mismatch of on-site energies. The main idea is to switch 
it on by the remaining gadget in our toolbox, a dynamical 
spin-vibron coupling (3) for the ic-ion. For spin-independent 
drivings Aco^ = 0, the periodic modulation of the on-site ener- 
gies assists the tunnelings by absorbing energy from the driv- 
ing to overcome the on-site gradient [14]. In this work, we 
propose to exploit the genuine spin-dependence of the driv- 
ings to build a single-spin heat switch and a current probe. 

For the single-spin heat switch, the parameters of the spin- 



vibron coupling (3) are Vk: — 2 
which lead to il^R(») = -i[//f^T 
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ia<PK 'a>PK 

Here, the tunneling becomes a spin operator through 
the argument of the first Bessel's function jf^^J (ai. ) = 

^ PkIo^ Pk' 

JpKia^i (C^(l + (^p,))^ where ^k = A(0+ /2v„ [22]. We note 
*^t 4tJMJ\ Ipk) = 0, while j;tliaj,J\ tp,) ^ 0. There- 



fore, by controlUng the ic-spin state via microwave ;r-pulses, it 
is possible to switch on/off the heat current through the TQD. 
Note that different switches have also been exploited in [41] 
to control the entanglement dynamics in harmonic chains. 

For a Ramsey probe of the vibron current, we need a 
minimally -perturbing mapping of the current onto the fc-spin. 
This requires a bichromatic spin-vibron coupling (3) with spe- 
cific parameters [42]. The first frequency assists the heat flow 

Hl^l = UJf5A-p.+^---)^ €1 - -i24..3i(;r), 

such that the dressed tunneling becomes purely imaginary. 
This is crucial to devise the probe, since the second frequency 
leads to Hi, = jhlj^^^p^, where /p, = i(/;^_ +rfpj. In 
the limit ^^.2 — > 0, A/ ~ 4^k:,2/^, we get a Ramsey probe (8) 
for the current mean value (/!''') and fluctuations 5,vib,vib(0). 

pk: 'pk'pk 

The ability to measure these fluctuations is essential for the 
comparison between fermionic and bosonic currents through 
the so-caUed Fano factor ^ = 5,vib,vib {0)/2{rJ^) . For flie heat 

Pk Pk "'*■ 

current trough the TQD, we expect strong super-Poissonian 
fluctuations, which are most pronounced for symmetric cou- 
pling to the thermal leads and increase linearly with fii in the 
regime n^ ^ n^ [12, 43]. Unlike the sub-Poissonian fluctu- 
ations in electrical cuiTents, the super-Poissonian fluctuations 
in heat currents have not been observed yet. 

Conclusions.- In summary, we have outlined the imple- 
mentation of an ion-trap toolbox for quantum heat transport, 
which provides (i) analogues of thermal reservoirs, quantum 
dots and wires; (ii) a way of engineering on-site disorder and 
dephasing, and (Hi) noninvasive probes for both vibron occu- 
pations and currents. We emphasize that the incorporation of 
these functionalities significantly extends the possible range 
of experiments on heat transport. Laser-cooling the edge ions 
to either coherent or squeezed vibron states [44], may consti- 
tute other valuable supplements to our toolbox. We expect, in 
addition, interesting effects in the presence of non-linearities, 
such as the interplay with Mott insulators [21], thermal rec- 
tification [45], and structural phase transitions [46]. In a 



non-equilibrium transport version of the spin-Peierls instabil- 
ity [47], it will also be interesting to explore if a finite heat 



current can be correlated to the structural change. 
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I. SUPPLEMENTAL MATERIAL: 
CONTROLLING AND MEASURING THE QUANTUM 
TRANSPORT OF HEAT IN TRAPPED-ION CRYSTALS 

In this supplemental material (SM), we present a detailed 
derivation of the expressions used in the main text. Besides, 
we test its validity by comparing the analytical expressions to 
numerical results for ion-trap setups with realistic parameters. 
Therefore, this SM will also be useful to guide an experimen- 
tal realization of the ideas discussed in our work. The contents 
of this SM are organized as follows: 
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Appendix A: Trapped-ion toolbox for quantum transport 

The objective of this section is to present a detailed deriva- 
tion, supported by numerical simulations, of the three gadgets 
in our toolbox: the vibronic tight-binding model (1), the con- 
trolled dissipation (2), and the spin-vibron coupling (3). 



1. Tight-binding model for the vibrons 

Let us start by describing how the Hamiltonian ( 1 ) for the 
transverse vibrons in a linear ion chain is derived. The dynam- 
ics of an ensemble ofN = Na +Nt:+Nk ions of three different 
species/isotopes is described by the following Hamiltonian 



fi=LL 



Inin 



1 



-ma(Oayr: 






L E 



1 



(Al) 
where v € {xjjjz} labels the different directions of motion, 
the latin indexes /,y' G {1 • • A^} label the ions, and the greek 
sub-indexes a,/3 € {cr, t, k] specify the particular type of ion. 
Here, (Hay are the trap frequencies for each ion species along 
the different axes, and e\ — e^ /^Tteo is expressed in terms of 
the electron charge and the vacuum permittivity. 

When the trap frequencies fulfill (Oa.z ^ (iia.xi(i>a.y, the 
equilibrium of the trapping and Coulomb forces leads to 
a Wigner-type crystal at low temperatures, the so-called 
Coulomb crystal. The crystalline equilibrium positions lie 
along the z-axis r? ~ z? e,, and are given by the solution of 
the following system of algebraic equations 
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Here we have introduced the relative masses jj.a = nia/nia, 
and the anisotropy of the trap frequencies «■« = (0^ Jo}^ ,. In 
this expression, we write the positions in terms of a common 
unit of length z^ — la^i , where l^, = (gQ/mfjCO^ _)'''^. For 
linear Paul traps, these equilibrium positions lead to an inho- 
mogeneous crystal, such that the ions in the center are more 
closely spaced than those at the boundaries [1]. By segment- 
ing the electrodes, it is possible to introduce a quartic com- 
ponent in the axial trapping potential, such that the crystal 
becomes more homogeneous [2]. Moreover, with the advent 
of the so-called micro-fabricated surface traps, it is possible 
to design the distribution of radio-frequency electrodes, such 
that the ions are held in any desired lattice geometry above the 
trap surface [3]. Therefore, we will also investigate homoge- 
neous linear ion crystal for quantum heat transport. 

As customary, a Taylor expansion to second order in the 
small vibrations around these equilibrium position lead to a 
quadratic model, the so-called harmonic crystal [4]. In the 
particular case of the ion chain, the vibrations along each di- 
rection r,gv = '"f V + ^'"lav get decoupled [1]. In particular, 
for the transversal direction v = x, after setting (Oa^x = Ca 
hereafter, the Hamiltonian becomes 
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where we have defined the following Coulomb couplings 
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Table I. Vibrational parameters for each ion species 
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The coupling between the vibrations of distant ions can also 
be understood as the result of a dipole-dipole interaction Vdd °^ 
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between the effective dipoles di^ = eSr,^ 



induced by the vibrating charges [5]. By quantizing the vibra- 
tions via the creation-annihilation operators 
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we get the announced quadratic model for lattice vibrons 
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Figure 3. Vibronic quantum dynamics: Exchange of a vibra- 
tional quantum (i.e. vibron) between two distant CT =^^Mg+, and 
T = Mg+ ions confined by a linear Paul trap with an equilibrium 
distance of lO^m. The solid lines represent the vibronic numbers 
given by the original Coulomb Hamiltonian (Al) (red: (ni^) for 
Mg+, blue: {112^) for Mg+), whereas the symbols stand for the 
vibronic numbers given the effective tight-binding model (A8) (red 
squares: (ni^) for ^^Mg+, blue circles: («2i) for ^"^Mg^). 



where we have defined the tunneling strengths 
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The final ingredient to obtain the desired tight-binding 
model (1) is to neglect the terms in Eq. (A6) that do not con- 
serve the number of vibrons, which is justified when the trap 
frequencies are much stronger than the tunneling strengths [6], 
namely y,^y„ <C (Oa + COp- Using a rotating wave approxima- 
tion (RWA), the Hamiltonian (A6) can be expressed as the 
sum of the vibronic on-site energy term H^o, and the vibron 
tunneling Hyt. This gives rise to the vibronic tight-binding 
model //tb = ^vo + Hyt of Eq. ( 1 ) in the main text, namely 
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As demonstrated in recent experiments [5, 7, 8], this model 
allows for a controlled tunneling of vibrons between different 
ions, which can be used for quantum-information processing. 
We would like to note that for the longitudinal vibrations of 
the ions confined in linear Paul traps, it is not possible to ne- 
glect the terms that do not conserve the number of vibrons. 
In this respect, the analogy to electron transport via the tight- 
binding model is compromised, although the system would 
still allow for the study of quantum heat transport. 

In order to show that the approximations leading to the vi- 
bronic tight-binding model are satisfied under realistic con- 
ditions, we perform a numerical comparison of the dynam- 
ics under Hamiltonians (Al) and (A8). The typical orders of 
magnitude for the vibronic parameters are summarized in Ta- 
ble I. Let us consider a particular example of a two-ion crys- 
tal with the species a =^^Mg+, and T =^'*Mg+. The trap 
frequencies are {(Oax,(»ay,(»a7)/^T^ = (5,5,0.5)MHz, which 
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lO/xm, and to 

a vibron tunneling strength of 7i^2<,/2?r « 12kHz. We con- 
sider an initial pure state with a single vibronic excitation 
along the x-axis of the a =^^Mg+ ion. As a consequence of 
the tunneling, this vibron should be periodically interchanged 
with the neighboring x =^'*Mg+ ion. In Fig. 3, we observe 
the predicted periodic vibron tunneling. The filled circles 
and squares are the predictions of the approximate Hamilto- 
nian (A8), whereas the solid lines correspond to the original 
Coulomb Hamiltonian (Al). To obtain the dynamics numer- 
ically, we truncate the vibron Hilbert space to Wmax = 2, and 
consider the three vibrational axis (i.e. 6 vibronic modes) with 
Coulomb non-linearities taken up to 8-th order (e.g. a? ). This 
simulation shows that the approximations leading to the tight- 
binding model are very accurate for realistic parameters. 



2. Atomic degrees of freedom 

Once the vibronic Hamiltonian (A8) has been derived, let 
us turn our attention to the atomic degrees of freedom of the 
ion chain. The different ion species can be divided into two 
groups, depending on wether we want to exploit the coherent 
C or incoherent '2 (i.e. dissipative) dynamics of the ions. In 
particular, we select the ion species T G 3, and d, fc G £. 

For the T-ions, we select a dipole-allowed transition 
|g£^) ^ \e^^) [Fig. 4(a)], whereas for the {ffjic} ions, the 
dipole-allowed transitions from a two-dimensional ground- 
state manifold {|g,-^ 1), |gia,2)} "^ k,a) define a so-called A- 
scheme [Fig. 4(b)]. We will now describe the atomic dy- 
namics for any particular species. In the absence of any 
laser radiation, the master equation for the atomic density ma- 
trix, p = — i[//s,p] + ^(p). This contains a Hamiltonian part 
//s = H^ +H^ +H^, where each ionic species a — {t, (7, k} 



evolves under 
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and we have introduced the different atomic energies er\, 
and the index n^ = 1, na,nK G {I72} which depends on 
the dimension of the ground-state manifold. Due to the 
coupling with the electromagnetic environment, the dipole- 
allowed transitions are also characterized by a spontaneous 
decay rate Fa- Considering the recoil effects of the emitted 
photons [9], the dissipation is described by a super-operator 
^(P) = La ^a{p), where we have introduced 



-S7 n S; „ • +H.C. 



(AlO) 



Here, we have defined the raising-lowering operators for each 
transition 5+ „_^ = k/a > (^'a ,«« I = iS7a.,J^^ ^"d integrated 
(summed) over all different directions (polarizations) of the 
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In this 



expression, k„^Uk is the wavevector of the emitted photon 
along the unit vector u^^, whose modulus is determined by 
energy conservation k„^ = l{£ea ^£ga„a^- Additionally, we 
have introduced the unit vectors of the atomic dipole operator 
Ua,«a' which depend on the angular-momentum rules for the 
dipole-allowed transition, and thus on the polarization of the 
emitted photon ^ e {0,±1} = {TTjff*} [9]. If we assume that 
a quantization magnetic field is applied along the trap axis (i.e. 
z-axis), we have 



Ua,na 



Cz, q^O, Ua,na 



1 



V2 



^TCx-iev), ^ = ±1. 



Let us note that the first term in the above dissipator con- 
tains the coupling of the atomic and vibrational degrees of 
freedom of the ions, which is a consequence of the recoil 
kicks associated to the photon-emission processes. This ex- 
pression can be simplified further when the vibrations are 
much smaller than the wavelength of the emitted light, namely 
^ia =*"? +^r,-^, suchthatfc„„|ui:-5r,„| ^ 1 (i.e. Lamb-Dicke 
limit). We Taylor expand the dissipator (AlO), and make use 
of the following integrals for dipole-allowed transitions 



d^al^iPa, /d^aCv-U^t^O, 



(A12) 



where v,/i G {x,y,z} account for the vibrational directions 
along the unit vectors ev,e|,, 5v,ij is the Kronecker delta, and 
we have introduced the following constants 
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Figure 4. Atomic level scheme for the different ions: (a) Two- 
level scheme {k£j),|e^^)} for a dipole-allowed transition with de- 
cay rate F^ of the T-ions, which is driven by a laser in a standing- 
wave configuration L™, such that the Rabi frequency is iiL™ • The 
standing wave is red-detuned Alsw < from the atomic transition, 
such that we can use it for laser cooling, (b) Three-level A-scheme 
{ I .?! a 1 ) ' k'a -2 ) > kin ) } foi a dipole-allowed transition with decay rate 
Fa of the a = {<y, K:}-ions. We use a couple of laser beams in a 
traveling-wave configuration, such that their Rabi frequencies for 
each of the optical transitions are ft/^„^, where «« = 1,2 stands 
for the two possible ground-states, and /« = lc(,2c( stand for the 
two laser beams. In this case, the corresponding detunings A/^ „^ 
will be much larger than any other scale of the problem, such that 
we can preserve the coherence of the ground-state manifold, and 
use the two-photon transitions to manipulate the state of the ion in 
this ground-state manifold by tuning the effective laser frequency 
G^L™ = fUi„ — tt>2a of the two beams. 



Note that, in analogy with the Coulomb couplings (A3), the 
recoil events to second order (A 12) do not couple the vibra- 
tions along different directions either Therefore, we can focus 
on the transverse vibrations along the x-axis directly (A5). 

Using these integrals, we can rewrite the dissipator (AlO) 
as a sum of two terms i^a(p) = i^J^(p) + ^^(p), where 



(All) ^^{•)= y -^(S; , •5+„ -5,+ „ S; „ .+H.C. 



(A14) 
describes the spontaneous emission of a collection of atoms 
with mutual distances much larger than the wavelength of the 
emitted light. In addition, the recoil effects are contained in 



^^(.) = V lfa5r„ {{a] +a. )*{a] +«• ) 
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H.c, 



where the strength of these terms fa = TaT]^ aqx is smaller 
than the bare dissipation (A 14) since Tj„^ = k„^/^2ma(i>a <C 
1. According to this expression, the photon recoil leads to 
dissipative events where the number of vibrons is modified. 

Typical orders of magnitude of some of these parameters 
are listed in Table II. In order to be more precise, let us con- 
sider a particular mixed ion crystal with species a =:2^Mg+, 
K =^Be+, and T =2^Mg+. The internal states corresponding 
to the level structure in Fig. 4 can be expressed in terms of 
the hyperfine atomic levels \nLj,F,M), where « is the prin- 
cipal quantum number, L,J are the orbital and total elec- 



Table II. Atomic and laser parameters for each ion species 
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tronic angular momenta, and F.M are the total angular mo- 
mentum and its Zeeman component along a quantizing mag- 
netic field. The T =^"'^Mg+ ions have no nuclear spin, and 
thus no hyperfine structure. For the two levels in Fig. 4(a), 
we choose je^j) = \iP\12) , l^f^) = |25i/2), such that the tran- 
sition frequency is C0q/2k = (gg^ — eg^)/2K ~ lO^THz, and 
the natural linewidth T-^/ln = 41.4MHz. Conversely, the 
O =^^Mg+ and K =^Be+ ions display a hyperfine struc- 
ture, which allow us to select two states from the hyperfine 
ground-state manifold and a single excited state, in order to 
form the desired A-scheme of Fig. 4(b). For o =^^Mg+, we 
take \gi^,i) = \3Si/2,2,2), \gi^2) = |35i/2,3,3), and the ex- 
cited state in the |e/^) = I3P3/2) manifold. The correspond- 
ing transition frequency between the ground-states lies in the 
microwave regime co" — (e^^i — £^^^)/2n ~ 1.8GHz, and 
there is a negligible decay rate (i.e. T(j/2k « 10^'^Hz). 
Therefore, all the spontaneous emission occurs via transi- 
tions to the excited state, which has a natural linewidth of 
Ta/2n = 41.4MHz. Finally, for K =^Be+, the ground-states 
would be |^,-^_i) = |25i/2,l,-l),|g,'^,2> = |25i/2,2,-2) with 
a transition frequency (Oq = (eg^, — eg^^)/2n = 1.25 GHz, 
and also a negligible linewidth. In this case, the excited state is 
in the manifold |e, J = |2Pi/2), such that r^/2;r = 19.4 MHz. 
Once these realistic parameters have been introduced, we 
can test the validity of the effective dissipators in Eqs. (A 14) 
and (A15). In Fig. 5, we compare the dynamics of T =^^Mg+ 
predicted by these effective terms to that given by the com- 
plete master equation (A 1 0) including the photon recoil in all 
possible directions. To obtain the dynamics numerically, we 
truncate the vibron Hilbert space to Wniax = 2 for each direc- 
tion V =x,y,z (i.e. three vibron modes), and consider all non- 
linearities of the recoil jump operators in (A 10) by discretizing 
the solid angle of the photon kicks in a mesh of 400 points. 
As follows from the displayed agreement, the description in 
terms of the dissipators (A14)-(A15) is very accurate. These 
equations are the starting point of the following sections. 



3. Edge dissipation by Doppler cooling 

We move onto the derivation of the effective dissipation (2) 
of the T-vibrons. We will be interested in positioning these 
ions at the edges of the chain, such that they can act as reser- 
voirs for quantum transport [Fig. 1(a)]. To control the dissi- 
pation of the edge vibrons, a laser beam L™ is tuned close 
to the resonance of the dipole-allowed transition. The master 
equation for the total density matrix p is 




Figure 5. Recoil dissipation for a dipole-allowed transition: (a) 

Decay of the population of the excited state Pg^ for T = Mg+ and 
an initial state Pt(0) = |eT)(ex| ® |0)(0|, where |0) represents the 
vacuum of vibrons in all three directions of motion. The solid line 
represents Pe^ according to the full dissipator (AIO), while the circles 
represent the same population calculated from the effective descrip- 
tion (A14) and (A15). (b) Dynamics of the mean number of vibrons 
{n\^) along the x-axis or T = Mg+ and the same initial state. Due 
to the recoil of the spontaneously emitted photons, the vibron num- 
ber is slightly increased (i.e. recoil heating) until the excited state 
has completely decayed. The solid line represents {n\^) according to 
the full dissipator (AIO), while the circles represent the same vibron 
number calculated from the effective description (A 14) and (A 15). 



where the dissipators S!a — ^a + ^a are given by Eqs. (A14) 
and (A 15). Additionally, the laser-ion interaction is 

^L-r - -Li^rSl + d;Si^ ) ■ El™ in^ , , (A17) 



where we have introduced the laser electric field ELSw(r£^,f), 
and the atomic dipole d^. For reasons that will become clear 
later, we need cooling rates that are much stronger than the 
vibron tunnelings (A8). Therefore, the laser beam is ar- 
ranged in a standing-wave configuration [10], Ei^swir^^jt) — 
eT:ET:Cos{cO]j'«t)cos{k^ -Tc^), where C-c^tjOl™ are the po- 
larization, amplitude, and frequency of the laser, and Rlsw is 
the laser wavevector directed along the .jc-axis (i.e. direction 
of the vibrons). Besides, we consider that the axis of the ion- 
chain lies at the node of the standing wave. 

We introduce the Rabi frequency Q.]j-« — —E-cdt ■ Ex, 
the Lamb-Dicke parameter tJlsw = ^l™/\/2?MtOt, and the 
detuning Al.™ = oJl™ - (e^, - Eg J [Fig. 4(a)]. When 
liiLrM^U"! < «o = i^e, - %), and r]u^ < 1, we can ap- 
proximate the laser-ion coupling (A17) by the expression 

^L|w-E^£.K+<), F,, = -ii2L.™rj^w5+e-''^™'+H.c. 

(A18) 
Since we are working at the node of the standing wave, let us 
note that the component of the laser-ion interaction that would 
drive the carrier is exactly cancelled. Therefore, the only fun- 
damental constraint over the standing-wave Rabi frequency 
will be |i2Lsw I <c Oq, still allowing for high driving strengths. 
To derive the effective dissipation of the T-vibrons, which 
was introduced in Eq. (2) of the main text, the crucial point is 
to appreciate the clear separation of time-scales 



Viajp I > I i^u" ^Lr I , irrT?^.™ I « r„ 



(A19) 



-i[//,+//tb+//L™,p] 



(A16) 



which implies that the spontaneous decay of the atomic states 
of the T-ions is much faster than the dynamics introduced by 
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any other term. This allows us to partition the Liouvillian as 
follows ^ — ^0 + ^1, where the two terms are 






^Lr-P]+^T(p) + ^a(p) + ^^(p), 

(A20) 



and where the "tildes" refer to the interaction picture with re- 
spect to //o = Hs +//vo in Eqs. (A8) and (A9). We can elimi- 
nate the fast degrees of freedom of the T atomic states by pro- 
jecting onto the steady-state of ^o(Ms^s) ~ 0- This can be ac- 
complished by projector-operator techniques [11], which lead 
to a second-order master equation 



(A21) 
where ^, and ^ — \ — !^ are the projectors of interest, 
which correspond to l^x\*] = Ms^s ® Tr^^atl*} in this case. 
Since the ion chain lies at the node of the standing wave, the 
atomic steady state is [il^ = ®i^ 1^^^) (^^'t I- After some algebra, 
Eq. ( A2 1 ) leads to an effective master equation for the reduced 
density matrix after tracing over the atomic degrees of free- 
dom of the T-species jU = TrT,at{p}- Back in the Schrodinger 
picture, this master equation reads as follows 

(A22) 
Here, we have introduced a dissipation super-operator that 
only acts on the vibrons of the ion chain, namely 



^(•) 



= A 
-A 



+ iJ 






•'/:," 



^ix' 



' LI // Li. I 



.)+H.c. 



(A23) 



As announced previously, the effective dissipation (A23) acts 
directly on the T-vibrons. The heating-cooling coefficients 
can be expressed in terms of the power spectrum of the laser- 
induced couplings F^^, namely 



(A24) 



F,^{(0)= / df(/>,(0/iaO))s 



In particular, the cooling depends on the power spectrum at 
positive frequencies A^ ~ Sf^ ^p; (+fi'('T)' while the heating 
is quantified by the noise spectrum at negative frequencies 
A^ = Sf( .f, i^^tx)- Such coefficients coincide with those 
of a single trapped ion [10], and it is the possibility of con- 
trolling experimentally the frequency asymmetry of the power 
spectrum which allows for an effective laser cooling of the vi- 
brational modes, i.e. Sfp ,F( {+(0e^) > Sf, .Ft (~fi^/'t)- By the 
quantum regression theorem [12], the coefficients can be ex- 
pressed as follows 






(jH^wTJliw) 

+ i(-AL|w±a)^J 



(A25) 



Note that in the expression above (A23), we have neglected 



rapidly-rotating terms (e.g. aj ), 



valid when o^^ ^ | A 

i^[A,0i,02](. 



= A(Oi.02-020i. 



an approximation that is 
v)^/2rT. Finally, by using 




Figure 6. Damping of the vibrons by laser cooling: Decay of the 

average number of vibrons («ij) for a single laser-cooled T =^'*Mg+ 
ion. The solid line corresponds to the predictions of the original 
master equation (A27), whereas the circles are given by the effec- 
tive dissipation (A23). We also display in a dashed straight line, the 
steady-state vibron number predicted by Eq. (A29). 



the dissipator (A23) corresponds to Eq. (2) in the main text. 

In Fig. 6, we compare the dynamics given by the effective 
edge dissipator (A23) with that given by the original master 
equation (A 16) restricted to a single T =^"'^Mg+ ion, namely 



p = -i[//,VHvo +//l™,p] + i^T(p). 



(A27) 



We consider an initial state Pt(0) — \eT:){eT:\ (8)p*, where p* 
is a thermal state for the T-vibrons with an average vibron 
number of n i j. = 1.5. In addition to the and atomic parameters 
for the T =^^Mg+ ions introduced above, we consider a trap 
frequency of (Ox/ln = lOMHz, and a standing-wave laser that 
is red-detuned Alsw = — jF^, such that its Rabi frequency is 
Hlsw = 0.1|Alsw|. In the figure, the solid line corresponds 
to the average number of vibrons (nir) given by the original 
master equation (A27), whereas the circles correspond to the 
effective description of the edge dissipator (A23). In this case, 
we need to truncate the vibron Hilbert space to Wmax = 15 to 
account for the thermal effects accurately. As follows from 
the agreement shown, the effective description (A23) in terms 
of heating/cooling processes at the edges is very accurate. 

Let us write down explicit expressions for the cooling rate 
Ye^ = ReKA^r )* — A^ }, and the mean number of vibrons in 
the steady state n^^ = Re{A^ }/7Cx' which read as follows 



7tx 



and 



i(nL|wT7Lsw)2rT,a)£^(-A^> 



((^r. 



" (Al™ 



-«£j2)((ir. 



' + (Al|w 



n<?T 






(A28) 



(A29) 



H.c. 



(A26) 



Note that the cooling rate (A28) changes sign depending on 
the laser detuning. In Fig. 7(a), this rate is represented as 
a function of the detuning in the so-called Doppler-cooling 
regime F^- ^ 0)^ of T —^'^ Mg+. For red detunings Alsw < 0, 
we get an effective cooling of the T-vibrons, whereas heating 
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Figure 7. Doppler cooling parameters: (a) Effective laser cooling 
strength yi^ for T = Mg+ (A28) as a function of the standing-wave 
detuning An* for a Rabi frequency i^L™ = sTt- For red detunings, 
we obtain cooling rates that can be as high as tens of kHz. (b) (left 
panel) Quadratic increase of the cooling rate as a function of the 
Rabi frequency, (right panel) Ratio of the nearest-neighbor vibron 
tunneling and the effective cooling rate J^^ (;+i)„/>/r as a function 
of the standing-wave Rabi frequency, (c) Steady-state mean vibron 
number ng^ (A29) as a function of the standing-wave detuning. 



is obtained for blue detunings Al^w > 0. Another important 
property is that the cooling rate (A28) increases quadratically 
with the laser Rabi frequency. In the left panel of Fig. 7(b), 
we see how the red-detuning cooling rate increases with the 
Rabi frequency without saturation, which shall allow us to at- 
tain regimes where the cooling is much stronger than the vi- 
bron tunnelings y^^ 3> |-/(a;g|, and the T-ions act as vibronic 
reservoirs for the heat transport along the ion chain. In the 
right panel of Fig. 7(b), we represent the ratio of the vibron 
tunneling and cooling rate as a function of the standing-wave 
strength. As announced previously, the cooling dominates 
over the tunneling for strong-enough standing waves. Finally, 
let us also note that the mean number of vibrons in the steady 
state (A29) is independent of the Rabi frequency. There- 



fore, increasing the laser power such that the desired regime 
Ye^ ^ \JiaJB I ^^ attained, does not limit the tunability over the 
vibronic reservoirs (see Fig. 7(c)), a property that will be im- 
portant to study the consequences of heat transport. 



4. Tailoring the spin-vibron coupling 

The final ingredient of our toolbox is the coherent spin- 
vibron coupling (3) for the ions a, /« G 2I- In particular, these 
ions will be positioned at the the bulk of the chain, such that 
the spin-vibron coupling can be used to control and measure 
the quantum heat transport [Fig. 1(a)]. In contrast to the 
engineered dissipation (A23), we use laser beams L'^ in a 
traveling-wave configuration that are far from the resonance 
of the corresponding dipole-allowed transitions [Fig. 4(b)]. In 
this section, we will use the indices a,j5 = {cT, K"} to refer 
only to the bulk ions. The corresponding laser-ion interac- 
tions //^tw are introduced in (A22) as follows 

A = -i[E(//s"+^L.w)+//tb,M]+E^^(M)+Li^a(M), 



(A30) 



where the laser-ion Hamiltonians for the A-scheme are 



//, 



(da,«a'5,„,„„ 



■da,„„'S;„,„„)-ELtw(r,„,f) 



(A31) 
Here, the electric field for each laser-beam arrangement is 
El'wK,?) = !/„£/„£/„ cos(k,„ -r/^ -0),„r), which consists 
of two la = la, 2a travelling waves with polarization, ampli- 
tude, and frequency C/^ , Ei^ , CO/^ , and a wavevector k/^ . 

To derive the spin-vibron couplings (3) from (A30), let us 
define the detunings and Rabi frequencies for each transition 



A/„„„ = «/„ - (ee„ - £; 



'Sajia I 



^lana = ^^l^d, 



a,na 



(A32) 

as depicted in Fig. 4(b). In the regime of weak excitations 
l^iasj I ^ ^lasi ' the auxiliary state is seldom populated, and 
we can obtain an effective description by projecting onto the 
atomic ground-state manifold. This is performed by using 
Eq. (A21), where the Liouvillian is now partitioned into 

^0(A)=E^a(A), 
a 

a f'j a 

(A33) 

and the "tildes" now refer to the interaction picture with re- 
spect to Ho — H^ + H^ +//vo- In this case, the projection 
operator ^^^ = L«„«^ \8o,n„,8K.n^){8a.na,8K.nJ restricts the 
dynamics to the ground-state manifold. After some algebra, 
and making use of the adjoint master equation for atomic op- 
erators [12], equation (A21) leads to an effective master equa- 
tion 



dA 
dT 



i[//tb,A] +E^v'(A)+A^a(A)+A^;c(A), (A34) 
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where p. — {^^ukP-} is the projected density matrix. The new 
super-operator in the LiouviUian can be expressed as follows 

A^a{')^Y.L L ^,X'"(r,„,OI^<a,«a)fe„yJ-+H.c., 

ia I a -I a "a, "'a 

(A35) 

where we have introduced the following couplings 

^"«n'a^_ f) _ ^^lana^l'gn'a -i((k,„ -k,, )t,„ -(A,„„„ -A,, , )r) 

(A36) 
Note that the form of these couplings (A36) tells us that the 
dynamics is due to two-photon processes that connect the 
ground-states via the excited state (see Fig. 4(b)). The real 
part of the coupling gives rise to an spontaneous decay due 
to the finite lifetime of the excited state, whereas its imagi- 
nary part leads to coherent dynamics within the ground-state 
manifold. As announced previously, in contrast to the en- 
gineered dissipation (A23), we use far-detuned laser beams 
Fa ^ |A/a»a I- In this regime, the couplings become 

(^"an'a(^ .n ^ ■ ^Igna^lgn'a -i((k,„ -k,, )-r,„ -(A,„„„ -A,, „, >) 
\l'a ^''-^')~' 4A,,„, ^ ' 

(A37) 
such that the spontaneous decay is minimized, and the coher- 
ent Hamiltonian part prevails leading to 



A^a(A) = -i[A//a,A]. 



(A38) 



Let us remark that the effective decay rates within the ground- 
state manifold scale as FaH^ia.nal / ^la.na)^ ^ T"' ^^'^ can be 
thus neglected for large-enough detunings. This is precisely 
the regime where F^ ;:^ Fcr, F^- considered in this work. 

By tuning the two-photon frequencies (O^jw = Oi^ — a>2„, 
such that jo^tw | <C {Sg^ , — ^nai) ~ ^o*' ^^^ lasers do not pro- 
vide enough energy to drive a two-photon Raman transition. 
Instead, they lead to the following two-beam Stark shift 

AHa=Y.Y. KT'^'"''^'^^Sia,na){8ia,J+ii.C.. (A39) 
"ataJ'a 

There are two kinds of terms in this expression AHa = AH" + 
AH". The two-photon processes whereby a photon is ab- 
sorbed from and re-emitted into the same laser beam (i.e. 
la — I'a) contribute with a standard ac-Stark shift 

AH"^^Y^Y^-Aea,„a\gia,na){8ia,na\> ^^a^a =2- 7\""" ' 
ia "a ^ la '^'^'aiia 

(A40) 

Additionally, when the photon is absorbed from and re- 
emitted into different beams (i.e. /« ^ l'^), the corresponding 
Stark shift leads to a more interesting term that couples the 
internal and vibrational degrees of freedom 



la "a 



where we have introduced the crossed-beam Rabi frequencies 



i2, 



■a ."a 2A- 



2ana 



(A41) 
gncies 

(A42) 



and wavevectors k^tw = ki^ ~k2„, and used (o" ^ |A/„„£,|. 
This crossed-beam Stark shift (A41) can lead to a variety of 
spin-vibron couplings. We discuss now how to produce the 
desired spin-vibron couplings (3). 

Let us introduce a notation where the two states in the 
ground-state manifold are referred to as spin states |§,„,i) = 
lt;„),l^<„,2) = |;,-J. We substitute r,„ = r^ + 5r,„ in 
Eq. (A41), such that the small vibrations are expressed in 
terms of the creation-annihilation operators (A5). By expand- 
ing in Taylor series of the Lamb-Dicke parameter 



^L'» - 



«1, 



(A43) 



" ^/IniaCOa 
and setting \Q.jjv, rj^jw | ^ Oa, we obtain 

A//« « U^La,,, (1 - 7]5.<a,Je""C'|,,J(.,J +H.C., 

(A44) 
where «,■„ — t, , U ■ The first term of this expression con- 
tributes with a periodic modulation of the Stark shift (A40) 

Ae„,„„ ^ Aea.na + I^L'»,„a I cos(coL.wf - 0L'„"), (A45) 

where we have introduced the phase of the Rabi frequencies 

ii, tw „ = |i2, tw „ le '-'a . The second term leads to 

La -ria I La "a ' 

^s"(0 = -EE I^Liy,.v,„ l^5„»'^°s(C0Ltwf-(/)Ltw)|i;„)(i,-JaJ^a,„ 
ia -^la 

(A46) 
We are now ready to derive the expression (3) used throughout 
this work. Let us make the following definitions 



together with the frequency and phase of the lasers 



Va = %tw, 9a = 0LS 



(A48) 



Then, the crossed-beam Stark shift (A46) becomes exactly the 
desired spin-vibron coupling in Eq. (3) of the main text. Let 
us note that the above drivings in the spin-independent regime, 
AcOa = 0, were used in [ 1 8] to mimic the effects of an external 
gauge field in the dynamics of the vibrons. 

Let us now support numerically this derivation. Note that 
the dynamics of the A-system (A30) for a single G 
ion, which is given by 



_25 



Ms 



A = -i [//f + //l;» + //vo ,ix]+^a{fi), (A49) 

covers a wide range of time dependences IGHz-lO^THz. 
Since we apply two laser beams with different frequencies, 
it is not possible to find a suitable picture where the above 
LiouviUian becomes time-independent. As a consequence, 
the numerics with the ^^Mg+ parameters require prohibitively 
large times. To avoid this numerical limitation, we shall work 
with dimensionless units, ensuring that the constraints for the 
derivation of Eq. (A39) are met. We set the atomic energies 
'■° ^8<t2 ~ ^'^ga 1 — ^'^ea ^ 50, and the laser frequencies to 
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Figure 8. Spin-vibron coupling: (a) Dynamics of the spin co- 
herence of a single ion initially prepared in the state |*I'(j(0)) = 
(l^tT.l) + \ga.2))/^f^ in the A-scheme [Fig. 4(a)] (see text for the 
particular parameters). The lasers are set according to ki^,k2„ || Cy, 
such that there is no spin-vibron coupling, and thus the lasers only 
induce a periodic oscillation of the coherences. The oscillation of 
the coherences {d'D (solid line: dynamics given by the master equa- 
tion (A49), whereas the circles correspond to the effective crossed- 
beam ac-Stark shift (A50)) confirms the analytical prediction, (b, 
c) Spin dynamics for the same state, but with the lasers fulfilling 
(ki^ — k2^) II e.v, such that the effects of the spin-phonon coupling 
become appreciable. In (b), we consider an initial vibrational Fock 
state with ha = (red solid line describes the coherences {d\) given 
by (A50); red circles correspond to (A44)), which coincide with the 
periodic oscillations of part (a). We also display ntj = 10 (pink solid 
line describes the coherences (d'j') given by (A50); pink squares cor- 
respond to (A44)), where one observes a shift of the oscillation pe- 
riod due to the vibronic state. In (c), we consider an initial thermal 
state with ha = W (pink solid line describes the coherences (d"j*) 
given by (A50); pink squares correspond to (A44)). In addition to 
the frequency shift, damping of the coherences is caused by the fluc- 
tuations of the number of vibrons associated to a thermal state. 



©i^ = 30, and 0)2^ = ft)i„ + 10"^ (e^^ ^ - %<j2)' ^^^^ "^^at the 
condition (i\jw = COi^ — ©2^ ^ COq leading to Eq. (A39) is 
fulfilled. The Rabi frequencies are set to iii^ 1 ~ —Q.2^^1 = 
O-Olggj,, and iii^ 2 = ^2a.2 ~ 0.016^^., such that the constraint 
|n/„.. I ^ A/„5 is also fulfilled. The final, and more relevant 
part, is to set the decay rate consistent with the far-detuned 
constraint Fa ^ |A/^„^|. We set Fa = lO^^Ai^j.i, which 
would correspond to detunings in the range 10-lOOGHz for 
the ff =^^Mg+ ion. Let us note that for these parameters, 
the standard ac-Stark shift (A40) is compensated, such that 
the dynamics should be given by the crossed-beam ac-Stark 
shift (A41). For our choice of Rabi frequencies, it becomes 



AH° 



n I2i tw 0^1 e 



-i(k,twr,CT-(B, twO 



H.c. 



(A50) 



where ii^- = 1.25 • IQ-^e,,, dj; = It Jit J - \UJiUJ- 
In order to show that this is correct, let us analyze the 
problem sequentially by first aligning the laser wavevectors 



kicr , k2^ II Cy, such that there is no coupling to the vibrons un- 
der study (i.e. x-axis). Since Eq. (A50) corresponds to a G' 
periodic driving, it should induce a periodic oscillation in the 
coherences of the qubit. In Fig. 8(a), we represent the spin 
coherences 

{&lit))^Tr{ia+^e-^<' + a^y<'Mt)}, (A51) 

where a^ = ||j^)(|j^| = (cTj" )1^. The solid line corresponds 
to the complete Liouvillian (A49), whereas the circles are 
given by the crossed ac-Stark shift (A50). From the agree- 
ment displayed in this figure, we conclude that the "adiabatic 
elimination" of the excited state leading to Eq. (A50) is very 
accurate. This allows us to work directly in the ground-state 
manifold, where the numerical limitations disappear, and we 
can use again the real parameters for the the a =^^Mg+ ion. 

Once this has been shown, let us align the laser wavevec- 
tors such that ki^ = —ii2a II ^x- In this case, the crossed-beam 
Stark shift (A50) introduces a coupling between the spin and 
vibrational degrees of freedom that will affect the coherences. 
We want to assess numerically the validity of the leading spin- 
vibron coupling derived in Eq. (A46). Hence, we consider 
a slowly oscillating o\jv, — lO^^cOa laser arrangement with 
i^Ltw = Q.lcOa/i'Hijw), where we recall that the transverse trap 
frequency for a =^^Mg+ ion is (Oa/2n = 5 MHz, and the 
Lamb-Dicke parameter is Tj^tw w 0.15. We truncate the vi- 
bron Hilbert space to Wmax = 60. In Fig. 8(b), we represent 
the spin coherences. We consider two initial Fock states with 
Ha G {0, 10}. First, we observe that the effective spin-vibron 
coupling (A46) is an accurate description. Second, we see 
that the period of the coherence oscillations depends on the 
number of vibrons, a feature that will be crucial to use this 
coupling as a measurement device. Finally, in Fig. 8(c), we 
initialize the vibrons in a thermal state with n^ = 10. We ob- 
serve that, for thermal states, the intrinsic fluctuations in the 
number of vibrons lead to a decoherence of the spin states. 
This feature will be crucial for heat transport measurements. 

In this appendix of the supplemental material, we have pre- 
sented a detailed derivation of Eqs. (A8), (A23), and (A46), 
which form our toolbox for quantum transport: the driven 
dissipative spin-vibron model (4). Moreover, all the different 
approximations have been supported by numerical analysis, 
which show that our analytical description is very accurate for 
realistic trapped-ion parameters. In the remaining part of this 
appendix, we will give details of the application of this tool- 
box for quantum heat transport. To ease notation, we write 
A — i^oKp-) -^ l-i in Eq. (4), bearing in mind that the atomic 
levels of the T-ions have been traced out, and we have pro- 
jected onto the ground-state manifold of the {<7, k} bulk ions. 



Appendix B: Thermalization: vibron number and current 

The objective of this section is to present a detailed deriva- 
tion, supported by numerical simulations, of the effective dis- 
sipation of the bulk vibrons (5), which forms the basis to un- 
derstand the ballistic heat transport across an ion chain (6)-(7). 
Additionally, we describe how to introduce dephasing and dis- 
order in the ion chain, and how they affect the transport. 
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1. Effective dissipation of the bulk vlbrons 

Let us derive the effective thermalization of the bulk vi- 
brons (5) starting from the driven dissipative spin-vibron 
model (4). In Fig. 7, we showed that the Doppler cooling by a 
standing wave leads to cooling rates yi>^ = Re{(A^ )* — A^ } 
that can be much stronger than the vibron tunnelings 2y(^ ^ 
Ji^j . In this limit, there is again a separation of time-scales: 
the thermalization of the edge T-vibrons is much faster than 
any other term in the Liouvillian (4). This allows for the par- 
tition of the Liouvillian into two terms 






(Bl) 



where the "tildes" refer to the interaction picture with respect 
to Ho = i/f + H^ + //vo- Let us start by switching off the 
spin-vibron couplings. To integrate out the edge vibrons, we 
use the projection-operator techniques (A21) for a projector 
■^edge{»} = Ms^s® TrT.vib{»}- Here, jj.^^ is the steady state of 
the laser-cooled T-vibrons determined by the different laser 
cooling parameters applied to each edge of the chain jj.^^ — 



nf' (g) /i* . In particular, it corresponds the thermal states 






-{l+nc. 



\ntr){nir 



(B2) 



with different mean vibron numbers n^^ = Re{A^ }/7ft- ^^ 
discussed in the main text, as long as the laser-cooling is 
switched on, the edge ions remain in a vibrational thermal 
state that can be controlled by the laser parameters. These 
edge T-ions act as a reservoir of vibrons for the remaining bulk 
of the ion chain. Our objective now is to derive the effective 
Liouvillian for the bulk ions. 

After some algebraic manipulations, and making use of the 
quantum regression theorem, we obtain the two-time correla- 
tion functions of the vibrons 



K(.)4, (0)),, = 5,.{n,^ + l)e-(rr^-'^^v)^, 



(B3) 



where we have introduced Ji^ — Re{(A^, )* — A^ } > 0, 5^,. — 
— Im{(A^)* — A^} > 0, and 5(^(1 is the Kronecker delta. 
Using these expression, together with the projection-operator 
formula (A21), we arrive at a master equation that only in- 
volves the degrees of freedom of the bulk ions 

Abuik = -i[l^i/s" +^^^Aibulk] + Aif (Mbuik), (B4) 
a 

where ^buik = TrT.vibJM}, and/Zj'^""' is the vibron tight-binding 
model restricted to the bulk ion species a, j3 G {(J, fc}. In the 
expression above, we have introduced the super-operator 



(B5) 



which is expressed in terms of the couplings 



•"JP %. -i((«<a-5f,) -«£,)■ 



(B6) 



The imaginary part of the T-coefficients can be rewritten as 
a Hamiltonian term, which leads to a renormalization of the 
vibron tunnelings and the on-site energies 

This leads to the renormalized tight-binding model 

^rtb = 1^ «;a<„a,a + 1^ E 47,3 «L«;,3 > (B8) 

introduced in Eq. (5) of the main text. In addition, the real 
part of the T-coefficients leads to a dissipative super-operator 



^buik(-) = EEfe/^ 



^ 'fl ; • fl; — fl; fl ; • ' 

Jp 'a 'a Jp I 



-A; ; {a- •«; —a; a- •) >+H.c., 

'ajp \ Jp 'a 'a Jp J \ ' 



(B9) 



where the dissipation rates are the following 



Kjp =ERe{<^,Kv, Kjp=L^^v^tjp}i"^r + i)- 

('t tx 

(BIO) 
Using the super-operator (B52), the above dissipator can be 
written as the bulk dissipator below Eq. (5) of the main text. 



2. Mesoscopic transport in ion chains 

The objective of this section is to provide numerical evi- 
dence supporting Eq. (5) for the effective thermalization of the 
bulk vibrons. Additionally, we will also check the accuracy of 
the predictions derived thereof, namely Eqs. (6) and (7) for 
the vibronic number and current through the ion chain 

Tl + Fr °- Fl + Fr 

(BID 

Following the philosophy of ''one, two, many", we first con- 
sider the smallest setup, a single-ion channel that will play 
the role of a thermal quantum dot (TQD). This will allow us 
to test the validity of Eq. (5), and the predicted number of 
bulk vibrons in the steady state. Then, we will move to a 
two-ion channel that will act as a double thermal quantum dot 
(DTQD). This will allow us to test the validity of the predicted 
current of bulk vibrons. Finally, we will explore the meso- 
scopic limit of a thermal quantum wire (TQW) formed by a 
longer ion chain, or a TQD connected to two thermal leads, 
where the leads are formed by a large number of ions. 

In order to test these predictions, we compare the dynam- 
ics given by the bulk master equation (5), with that given 
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by the master equation where dissipation only occurs at the 
edges (4). Since both Liouvillians are quadratic in creation- 
annihilation operators, it is possible to obtain a closed system 
of (A^ — 2)^ or A^^ differential equations for the two-point cor- 
relators Ci^j„ = (flj aj ), respectively. Both theories can be 
recast in a matrix differential equation 



a 1 



dC 
dT 



:i[J,C]-(WC + CW*)+I 



(B12) 



where the matrices J,W,K depend on the particular master 
equation, Eq. (4) or Eq. (5). In particular, for the edge dissi- 
pation (4), we get the following matrices 



'ajp 



((A, 



-K)^ia.eApA: «.^ e {a, k,t} 



<S=^-{K}^..iApA^ 



whereas for the effective bulk dissipation (5), we get 



(B13) 



jjbuik ^ 2Re{A+ ,. }. 



a,j3e{(j,K-} (B14) 



The possibility of expressing the dissipative dynamics as a 
closed set of differential equations (B 1 2) allows us to circum- 
vent the numerical Umitations, which would arise due to the 
large truncation of the vibron Hilbert space required for some 
of the simulations of the dissipative vibron model. 



a. Thermal Quantum dot: vibron number 

Let us consider the minimal scenario: the thermal quantum 
dot. In this case, the chain is composed of three ions 
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Figure 9. Thermalization of a thermal quantum dot: (a) Dynam- 
ics of the vibronic numbers (n,-^\ for a Mg+- Mg+- Mg+ ion 
chain (see the text for the particular parameters). The solid lines rep- 
resent the numerical solution of Eqs. (B12)-(B13), showing that the 
edge ions thermalize much faster (see also the inset for (ni \, ("Sj))- 
For the bulk ion, the numerical solution for the vibron number («2 ) 
given by Eqs. (B 12)-(B 14) is displayed with red circles, and shows a 
good agreement with the previous dynamics, (b) Steady state number 
of bulk vibrons (n2„ ) as a function of the reservoir couplings Fl /Fr . 
The solid yellow line corresponds to the prediction (B 1 8), the green 
crosses are obtained from the numerical integration of Eqs. (B12) 
and (B13), and the red circles correspond to the numerical solution 
of the effective dynamics in Eqs. (B12) and (B14). The grey area 
marks the region where the constraint lyi^ 2> Ji^jg is not fulfilled 
any longer, and thus deviations from the theory appear. 



such that the heat transport takes place along the minimal 
channel, namely a single-ion connecting the two T-reservoirs. 
According to Eq. (5), the bulk vibrons evolve under 

Abulk = -i[{COa + J2a2„)al^a2^, fibulk] + ^2^2,, (^bulk), 

(B16) 
where the dissipator contains the heating and cooling terms 



^2<.2<,(»)=At2„(fl2,»«2a- 



l/t l/t " 



-Ao -, (a^ •a\ — a\ Qr. •' 

ZaAa V ^a ^a ^a ^a ' 



H.c. 



(B17) 



Such a single-oscillator master equation can be solved exactly, 
and yields a non-equilibrium steady-state mean vibron num- 
ber of n2<, = Re{Aj^2„}/Re{A2^2<, - Aj^2„}- According to 
Eq. (BIO), this mean vibron number is the following 



«2„ 



rL«L + TrHr 



(B18) 



where we have introduced the mean vibron numbers «l = 
nit;«R = «3r foi" the left and right reservoirs. Here, we have 



also used the couplings Fl 

were introduced below Eq. (5) in the main text, namely 






and Fr = F-,'^ , , which 



-.i, _ 



2(7 ,2o 



= 27r72<,i,Pi;,(fi)2„ )-/£„: 



(B19) 



where Pfj.(e) is the Lorentzian density of states for the laser- 
cooled ions. Let us note that, considering the symmetry of the 
vibron tunnelings /2nit — ■lia'^x^ ^^ ^'^'^ '^^at the steady-state 
vibron number of the TQD does not depend on the Coulomb 
couplings, and thus on its distance to the reservoir. 

Let us now consider the realistic parameters for a x-O-x 
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Ms 



In addition to the 



chain, where T = Mg+ and a 
parameters introduced in previous sections, we consider the 
detunings Ai^ = — 0.6Ft, As^ = — O.SF^, and the Rabi fre- 
quencies Olsw = iiLsw = Ft for the laser cooling of the T- 

1t 3t 
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ions, where we recall that YT:/2n = 41.4MHz. With these 
parameters, the effective cooling rates (A28) of the T-ions 
would be JiJlK « 86kHz, and y^j2n « 106kHz. Addi- 
tionally, the mean number of vibrons for each reservoir (A29) 
would be h\^ = 1.65, and «3^ = 1.63. The trap frequencies 
are {a)ax,0)ay,0)az)/'^^ = (5,5,0.5)MHz, which lead to an 



inter-ion distance of zV 



: 9/im, and to a 



vibron tunneling strength of Jx^^^l^'^ = htial'^'^ ~ 30kHz. 
With these parameters, the constraint 2Y(\ ^ /,„;„ is fulfilled, 
and we expect that the T-ions thermalize very fast and act as 
an effective reservoir for the bulk d-ion. 

In Fig. 9(a), we confirm this behavior numerically. The 
blue solid lines represent the vibron number of the edge T- 
ions, namely (ni^), («3^), obtained by integrating Eqs. (BI2)- 
(B13). As displayed in the figure, the edge vibrons thermalize 
on a jUs-scale (see also the inset). The red solid line stands 
for the bulk vibron number {112^), which reaches the steady 
state on a longer millisecond-scale. The red circles represent 
the numerical solution of the effective bulk dissipation from 
Eqs. (B12)-(B14), whereas the red solid line corresponds to 
Eqs. (B12)-(B13). The agreement of these results shows that 
the effective bulk Liouvillian (B4) is a good description of 
the problem. Moreover, the red dashed line represents our 
prediction for the stationary bulk vibrons (B 1 8), which also 
displays a good agreement with the numerical results. Fi- 
nally, the blue dashed lines represent the laser-cooling vibron 
numbers «i^ , n^^, which perfectly match the edge steady state. 
This shows that the T-ions behave as a perfect vibronic reser- 
voir, whose properties are not modified by the bulk ion chain. 

In Fig. 9(b), we display the dependence of the steady-state 
number of bulk vibrons («2cr) as a function of the laser-cooling 
parameters at the edges of the chain. We keep the same pa- 
rameters as above, but let one of the laser-cooling Rabi fre- 
quencies vary in the range Q.]jv e [O.IF^, lOF^], while the 
Other one is fixed Q.U-* = F^. This does not modify the reser- 
voir vibron number, but only the bulk-reservoir couplings to 
the Fl,Fr. In the figure, we again observe a regime with a 
very good agreement of the theoretical prediction (B18) and 
the numerics from Eqs. (B12)-(B13) or Eqs. (B12)-(B14). It 
is also very instructive to analyze the shaded region, where 
deviations from the theory start arising. In this region, the 
condition lyc^ 3> Jiajg is no longer fulfilled, and one cannot 
treat the edge ions as an effective reservoir any more. This is 
the underlying reason for the displayed differences. 



h. Double thermal quantum dot: vibron current 

Let us now move into the double thermal quantum dot, 
which is formed by a two-oscillator channel connected to the 
two laser-cooled reservoirs 

T - C7 - (7 - T, (B20) 

In this case, we choose 

= — 0.6Ft, and the Rabi 

. = Ft for the laser cool- 

ing. The trap frequencies are set to {(0ax,(0ay,(0az)/'^'^ — 



where T =^'*Mg+ and a =^^Mg^ 
the detunings Ai^ 
frequencies Cli^s^ 
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Figure 10. Vibron number and current In a double thermal 
quantum dot: (a) Thermalization dynamics for the number of bulk 
vibrons («2„), (i3„) (see the text for the particular parameters). The 
solid lines correspond to the numerical solution of Eqs. (B12)-(B13), 
and the symbols to the numerical solution of Eqs. (B12)-(B14). 
The grey dashed lines represent the reservoir mean vibron numbers 
nL,/JR, while the red dashed line stands for the theoretical predic- 
tion for the bulk vibron number (B 1 1). (b) Steady-state bulk vibron 
number (n2„) = ("3„) 3S a function of the system-reservoir effective 
couplings rL,rR. The green crosses represent the numerical solu- 
tion of Eqs. (B12)-(B13), the red circles that of Eqs. (B12)-(B14), 
and the yellow solid line corresponds to the theoretical prediction 
in Eq. (BU). (c) Same as above, but displaying the steady-state 
vibron currents (/2''^)s.s according to Eqs. (B12)-(B13) (crosses), 
Eqs. (B 12)-(B 14) (circles), and the prediction (B 1 1) (solid line). 



(5, 5, 0.2) MHz. We consider an initial thermal state, where 
the two (T-ions have a different vibronic number n2g (0) — 2.4, 
and n3g (0) = 1 . Accordingly, we expect to observe a periodic 
exchange of vibrons between the bulk ions, which is addition- 
ally damped due to their contact with the reservoirs. 

In Fig. 10(a), we show the thermalization dynamics of such 
a two-oscillator channel. The solid lines correspond to the nu- 
merical solution of Eqs. (B12)-(B13) {{n2„) yellow, and (nj,^) 
green), whereas the symbols stand for the numerical solution 
of Eqs. (B12)-(B14) {{112^) yellow crosses, and (n3^) green 
crosses). The clear agreement supports once more the validity 
of our analytical derivations. More importantly, we also ob- 
serve how both vibronic numbers tend to the same equilibrium 
value given by Eq. (Bl 1) (red dashed line). 

Let us now address the validity of the predictions (Bll). 
In order to calculate the vibron current, let us note that the 
current operator can be defined through a continuity equation 
dn;„/df = /^^ -7^'? . By applying this to the tight-binding 
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Hamiltonian (A8), we get 



rvib 
rvib 



-iV y J*: fl/ fl; +H.C., 

P jp>'a 



(B21) 



Jp>'a 



where we have used Ji^j^ ~ JjRia- ^'^ '^he particular case of 
Eq. (A8), the tunnelings are real. However, we keep the above 
expression general since it will be useful in other sections be- 
low. In Figs. 10(b)-(c), we let one of the Rabi frequencies 
vary in the range Hl^w G [O.IFt, IOFt], which allows us to 
modify the ratio Fl/Fr. As shown in these figures, when the 
constraint lyi^ ^ JiaJR is fulfilled, there is an excellent agree- 
ment of both numerical solutions based on Eqs. (B12)-(B13) 
(green crosses), or Eqs. (B 12)-(B 14) (red circles), and the the- 
oretical predictions (Bll) (yellow lines). As expected, we can 
thus conclude that Fourier's law is violated in the TDQD. 



c. Thermal quantum wire: assessing Fourier's law 

Let us now consider a mesoscopic TQW with N — 20 ions, 
which would have a length of L « 0.21 mm for the trap fre- 
quencies (oJa.v, Way,tOaj;)/27r = (5,5,0.1)MHz. The configu- 
ration of ions species is 



T, 



(B22) 



where T =^'*Mg+ and a =^^Mg+, and we choose the detun- 
ings Aij = — 0.8Ft, A^^ ~ — 0.6Ft. We shall use this setup to 
test the validity of Fourier's law of thermal conduction. As 
discussed in the main text, this law predicts the onset of a lin- 
ear gradient in the number of carriers between the reservoirs 



FL =nL- 



(»r-«l) 



(B23) 



In Fig. 11 , we represent the number of vibrons in the steady 
state of the TWQ when the laser-cooling Rabi frequencies are 
set to i^Lsw = iij^sw — 1.4Ft. The yellow bars represent the 

It lOr 

number of bulk vibrons {nt^ ) ss obtained numerically from the 
steady-state solution of Eqs. (B12)-(B13), whereas the green 
bars stand for the numerical solution of Eqs. (B12)-(B14). 
These numerical simulations confirm the theoretical predic- 
tion (B 11 ) to a good degree of accuracy. It is also clear from 
this figure that the number of vibrons does not display a linear 
gradient, as predicted by Fourier's law, but is rather homoge- 
neous along the whole ion chain. 

As mentioned in the main text, the violation of Fourier's 
law is not a surprise, since this law applies to diffusive pro- 
cesses, whereas our vibron transport is ballistic. Let us now 
explore two possible mechanisms to introduce diffusive dy- 
namics in the problem. 

/) Noise-induced pure dephasing- A possible mechanism 
to introduce diffusion in the transport of heat is to consider an 
engineered noise leading to dephasing in the vibron tunnel- 
ing. This can be accomplished by injecting a noisy signal in 
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Figure 11. Violation of Fourier's law: Steady-state number of 
vibrons (n,-^) along the ion chain (see the text for the particular pa- 
rameters). The bulk of the TQW displays a homogeneous number of 
vibrons, in clear contrast to the prediction of Fourier's law. The yel- 
low (green) bars correspond to the numerical solution of Eqs. (B12)- 
(B13)(Eqs. (B12)-(B14)). 



the trap electrodes [13], leading to fluctuating trap frequencies 
that modify the on-site energies of the tight-binding Hamilto- 
nian 



^vo -^ Hyo + 5//vo(f) = ^(«,„ + 8cOi, 



(0)flL«', 



(B24) 



Here, we have considered that 5cOi^ (f ) is a zero-mean random 
Markov process that is stationary and Gaussian. Such pro- 
cess, usually known as the Ornstein-Uhlenbeck process [14], 
is typically characterized by a diffusion constant c, and a cor- 
relation time Tc, which we assume to be much shorter than 
the time-scales of interest f ^ Tc. Moreover, we introduce a 
coiTelation length ^c in order to model the extent of the noisy 
signal on the trap electrodes. The power spectrum of this noise 
is defined as 



(B25) 



Ss(o,^ .Scoj (o)) = Re <^ / df 5a),„ (f )5a);^ (0)e 







where the "bar" refers to the statistical average over the ran- 
dom process. In particular, the above three constants deter- 
mine completely the noise spectrum 
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(B26) 



where we have introduced the equilibrium positions of the 
ions given by Eq. (A2), and the rate F^ — ct^/2. 

By using a Born-Markov approximation to study the effects 
of the fluctuating trap frequencies, the master equation for the 
TQW becomes 



A==%dsv(Ai)- / dt'[8H,o{t)ASH,o{t-t'),fi]]. (Bll) 
Jo 

Using the above noise spectrum, the Liouvillian of the TQW 
gets the additional contribution of a pure-dephasing super- 
operator ^ddsv -> -S^ddsv + ^d, where 

a,pia,jp 

(B28) 
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such that the dephasing rate only depends on the zero- 
frequency component of noise spectrum. We also observe that 
^c controls the collective effects in the dephasing dynamics of 
the TQW: when ^c -^ 0, we obtain a purely local dephasing 
that introduces phase-breaking processes in the vibron trans- 
port, whereas for ^c ^ °°, the noise is purely global, such 
that the tunneling dynamics is not affected, and remains in 
the purely ballistic regime. 

This collective dephasing modifies the system of differ- 
ential equations (B12) for the two-point vibron correlators 
C,„,j„ = {a]^aj ), which becomes 
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where we have introduced the following matrix 



LL 2rdl-e 



a.piajjl 



\l^ 



l'a)(;/3l, (B30) 



where Wia)Yi =\ form an orthogonal basis of the A^- 
dimensional subspace of the two-point correlators. 

In Fig. 12, we compute the steady state solution of the 
above system of differential equations (B29). We consider the 
same experimental parameters as previously, and set the de- 
phasing rate to Fd = 10')^^. As can be observed in Fig. 12(a), 
in the limit of large correlation lengths ^c ^ |zir ^^n^V the vi- 
brons display the same homogeneous distribution that violates 
Fourier's law (i.e. balhstic regime). As the correlation length 
is decreased, a linear gradient starts to develop at the bulk 
of the chain (diffusive regime). It is interesting that we have 
a single parameter to control the ballistic -diffusive crossover. 
However, note that edge effects mask the linear gradient. 

We have found that these edge effects are particularly 
strong for a linear Paul trap, since the equilibrium positions 
in Eq. (A2) lead to an inhomogeneous crystal. By modifying 
the dc trapping potentials, or by considering micro-fabricated 
ion traps, it is possible to obtain an homogeneous ion crystal. 
In Fig. 1 2(b), we study numerically the distribution of vibrons 
in this regime, as the correlation length is decreased. Our re- 
sults show that edge effects are less pronounced in this case, 
and a perfect linear gradient, as predicted by Fourier's law, 
develops along the whole chain. 

//) Spin-assisted random disorder- In order to introduce an- 
other diffusive mechanism for the transport of vibrons, we ap- 
ply a spin-vibron coupling (3). By controlling the laser in- 
tensities, polarizations and detunings, we further impose that 
A(0^ = 0, which leads to a static spin-vibron coupling 



//f, = XiAco„ 



(B31) 



whose strength Ao,^ can be controlled at will, and made 
strong enough. The idea to mimic the effects of diagonal dis- 
order is to use the spin degrees of freedom as a gadget to build 
a Liouvillian with random on-site energies. Here, the random- 
ness is inherited from the quantum superposition principle in 
the spin degrees of freedom [15, 16]. 
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Figure 12. The dephasing route to Fourier's law: (a) Steady-state 
number of vibrons {n-,^ ) for an inhomogeneous ion chain in a linear 
Paul trap. As the correlation length of the dephasing noise ^^ de- 
creases (in units of the nearest-neighbor spacing at the center of the 
chain), keeping Fj = lOyA'^. we observe a inhomogeneous distribu- 
tion of vibrons across the chain. Far away from the edges and close 
to the bulk of the chain, the distribution displays a linear gradient. 
(b) Steady-state number of vibrons (n,-^) for an homogeneous chain 
in micro-fabricated ion trap array. As <^c decreases and Fj = 10)5Vj 
is fixed, we observe a perfect linear gradient across the chain. There- 
fore, edge effects are less pronounced in this homogeneous scenario. 



Let us consider an initial pure state for the d-spins of the 
TQW, namely Mspm(O) = l^oX^ol- Without loss of general- 
ity, it can be expressed as j'Po) = £{j. }C{^.. jKii^}), where 
{■Sial = {■S2<,,«3a:--- :«A'-ia} J^ ^ particular spin configura- 
tion for the bulk a-ions Si^ G {tiojiin}- The reduced density 
matrix of the vibrons evolves in time according to 

Mvib(0 = Tr,p,n{e^''^-«'^^»Vspin(0) ® Mvib(O)}, (B32) 

where we have rewritten the spin-vibron Liouvillian (4) mak- 
ing explicit reference to its dependence on the spin operators 

-^dsvdo'^J') =-^ddsv({(^|„,o-|„,--- ,(^N-\„})- Fromthisex- 
pression, the reduced density matrix evolves as 



Mv,b(0= L /^{..^le^-i-^^-JVviblO), 



(B33) 



which can be interpreted as an statistical average of the time- 
evolution under a stochastic Liouvillian. In particular, the Li- 
ouvillian ^ddsvd'^icr}) depends on the binary variables {.?,„}, 
which inherit their randomness from the quantum parallelism 
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of the initial spin state. In fact, the associated probability dis- 
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tribution for the binary random variable is p^^. i 
Therefore, we can formally write /Xvib(f) 



-{^..}l 



Mvib(0' where the 
'bar" refers to a statistical average over a random Liouvillian 

=%dsv -^ -^sdtb (Mvib ) = -i [Hsth , Mvib] + Y. ^v' (^vib ) • (B34) 



Here, i^y^ is the dissipator acting on the edge vibrons (2), 
whereas the stochastic tight-binding Hamiltonian is 



^stb = 



LEi^a] a- 



a,pia^Jl3 



T 



+H.C.). 



Here, the on-site energies of the bulk d-ions are binary ran- 
dom variables sampling £,-^ e {o),-^ — ^Ac0fj,C0i^ + jAco^}. 
In particular, if we consider the initial spin state |^o) = 
•^iffdt!- ) + li; ))/\/2, this diagonal disorder has a flat prob- 
ability distribution p{£i„) = j- Let us note that exactly the 
same formalism would apply to the effective Liouvillian for 
the bulk vibrons. In this case, however, the collective dissipa- 
tion rates would also become stochastic variables. 

In order to study the steady state for the vibrons thermal- 
izing under this disordered Liouvillian (B34), we can solve 
the system of differential equations for the two-point correla- 
tors (B 12) for each realization of the diagonal disorder 



dC 



{£,<.} 



dr 



^i[J'q.,.}]-(wq,,^}+q,,^}W* 



(B35) 



Then, we should average over the random variable according 
to the probability distribution p{£i„) = I 
one of the matrices becomes stochastic 



to the probability distribution p{£i„) = j. Due to the disorder. 



^fXi{>^kJ)-'^iaSiaJp+J. 



Jp '''•a.lp ■ 



(B36) 



After performing the statistical average C{t), we can recon- 
struct the vibron density of the disordered TQW. 

We consider the same setup as in Fig. 1 1 for the ordered 
TQW, namely a A^ = 20 ion chain. Moreover, we use the 
same parameters introduced there. For the spin-induced dis- 
order, we set Aco^ = IO/at^, which corresponds to a strong 
spin-vibron coupling. In Fig. 13, we represent the distribution 
of vibrons along the TQW in the steady-state. In this case, the 
predictions for both a homogeneous ion crystal (i.e. microtrap 
array), and an inhomogeneous one (i.e. linear Paul trap) coin- 
cide. As a consequence of the disorder-induced diffusion, the 
vibron layout is no longer homogeneous, but rather displays a 
linear gradient far way from the edges of the chain. 

Before closing this section, let us also comment on another 
interesting perspective for the TQW, namely the possibility of 
realizing noise-assisted quantum heat transport. As demon- 
strated in [17], the efficiency of transport in quantum networks 
including linear chains with disorder may be sometimes in- 
creased by the presence of local dephasing noise. In order to 
test this prediction in our current scenario, let us note first that 
the presence of disorder (B34) will partially inhibit the heat 
transport. By switching on the local dephasing (B28), the in- 
terference leading to Anderson localization, or transport bot- 
tlenecks due to energy mismatches between neighboring sites. 
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Figure 13. The disorder route to Fourier's law: Steady-state num- 
ber of vibrons (n,-^) for an inhomogeneous ion chain in a linear Paul 
trap (purple crosses), and a homogeneous ion chain in a microtrap 
array (green crosses). We observe an inhomogeneous distribution of 
vibrons across the chain. Far away from the edges, and close to the 
bulk of the chain, the distribution displays a linear gradient. 



can be overcome thanks to the presence of noise, thus assist- 
ing the transport of heat. This can be probed by the current 
measurement described in a section below. 



d. Effective thermal leads and single-spin heat switch 

Let us now consider a mesoscopic ion chain with A^ ions. 
The configuration of ions species is 



■T, 



(B37) 



where T ^^'^Mg^, a =^^Mg+, and K =^Be+. The JC-ion plays 
the role of the thermal quantum dot (TQD), the T-ions act ef- 
fective vibronic reservoirs, and the left/right chain of CJ-ions 
acts as leads that connects the TQD to the reservoirs. We will 
start by discussing the conditions under which the ff-ions can 
be interpreted as effective thermal leads. Then we will dis- 
cuss how to control the tunneling of vibrons across the fc-ion, 
which can be exploited to build a single-spin heat switch. 

/) Effective thermal leads- In order to devise the leads, we 
apply a strong and static spin-vibron coupling (3) to the a- 
spins (see Eq. (B31)), namely 



//c. 



,5^®a«<a(^4' 



In this case, we consider the strong-driving regime 



(B38) 



(B39) 



and the following initial state for the spins of the leads 

|ro> = |;<T---;a)|^J|t<T---ta>, (B40) 

where {(p^) is an arbitrary spin state of the fc-ion acting as the 
TQD. In this regime, the spin-vibron coupling (B31) provides 
a large and static shift of the vibron on-site energies 

COi^ -^ (bi^ = COi^ - \A(0„e{p^-ia) + \A(i)^e{ia-pK). 

(B41) 
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where we have introduced the Heaviside step function (x) = 
1, if jc > 0. Due to these shifts, the effective thermahzation of 
the bulk ions described in Sec. B 1 must be re-addressed. As- 
suming that the separation of time-scales is valid, whereby the 
laser-cooled T-ions reach the steady state much faster than any 
other process in the problem, we can derive a similar master 
equation (B4) for the bulk ions 

Abuik = -i[Y,H,"+ //,'?,"'^ Mbuik] + A^(Mbuik) , (B42) 
a 

where a ~ {<7,K"} are the bulk species. In this expression, 
the tight-binding model for the bulk vibrons must be modified 
to account for the on-site energy shifts (B4I), and reads as 
follows H^^^^ = H^^^^ + H^^'^, where 



H, 



bulk 



dj,„a,' fl; , H^^^^^Y, E ■^'ai/j«L«--+H-c. 



"'oc'^ia'^ia' "vt 



a-.P ia^jf, 



^iria"jfi 



(B43) 

Additionally, the second-order processes where vibrons tun- 
nel between the edge and bulk lead to a similar super-operator 



a,/3(aj'pi'i: 



+ «<T(fl,„»fl;/i-«>p«L» 



H.C., 

(B44) 



where the u -couplings (B6) are modified due to the shifts 

w„ iajp r,,-i((a),„-5,g-«,j- ^ ' 

By assuming that all the ion species have the same trap fre- 
quencies, and that the conditions (B39) are fulfilled, one can 



check that rr ,. «T 



(r 



If we now move to the interaction 



'ajp lajp 

picture with respect to Hq = H° + H^ + //^""', we obtain 

Abuik - -i [H^fHt) , Abuik] + Ai?(Abuik) , (B46) 

where we have introduced the vibron tunneling Hamiltonian 



^vrw-E E J, 



-.//3<%e^ ^ +H-C-: (B47) 



«:/3 'a^ip 

and the super-operator for the virtual edge-bulk tunnehng 



ijp 



H.c. 

(B48) 



Note that, when the conditions (B39) are satisfied, we get 

ferent halves of the chain. Therefore, a rotating wave approx- 
imation allows us to neglect all the tunneling processes that 
lead to the thermahzation between the two halves of the ion 



chain. This observation allows use to rewrite the partition the 
master equation (B46) into three terms 

Abuik = ^L (Abuik) + ^L^cR (Abuik) + ^R (Abuik) • (B49) 

Here, we have introduced the Liouvillian for each of the leads 

^l/r(») = -i[^L/R,»] + ^l/r(»), where 



Hl= E4...«l«/.e+'(»--«->+H.c., 

•o.j<r<PK 



(B50) 



'a,Ja>PK 



where we have used the renormalized tunnelings of Eq. (B7). 
Additionally, the dissipators are 

^L= E ^[A:,,,«L(o,«,-<,(0]+i^[Ar.,;v«,.(o.4w]' 
^R= E ^[^:..'<w,«;<.(0]+^[Ar,,,„,«,.(0,4w]> 

'c!,ia>PK 

(B51) 

where we have used the dissipative couplings in Eq. (BIO), 
the interaction picture operators fl;^(f) = a,j,e^"*''<j', and the 
generic super-operator 

@[K, Ox , O2] (.) = A(Oi . O2 - O2O, .) + H.c. (B52) 

The final part is the coupling of the leads to the TQD, which 
reads 

=^l^r(») = -i[HL^R(f),»] +^^(»), (B53) 

where we have introduced the Hamiltonian 

/?l^r(0- E 24,P.«L«p,e+i('^'--'^'^)'+H.c., (B54) 

and the dissipator due to the long-range tunneling between the 
reservoir and the TQD 



^- = ^[^P.,P.'4.'«/.J + ^[V,P.'«P.'4J' (B55) 

From the master equation (B49), we thus expect that the 
left/right half or the chain thermalizes individually with the 
left/right reservok, such that the vibron number in the steady 
state is 



\ni„ 



■ ni^e{pK-ia)+nRO{i(,~PK)- (B56) 



Thus, the two chains of d-ions serve as a thermal lead to con- 
nect the reservoirs to the thermal quantum dot (i.e. the fc-ion), 
and modify the local density of states seen by the TQD. 

In order to support this theoretical prediction, we integrate 
numerically the system of differential equations for the two- 
point correlators C,„j„ — [a] aj„), namely 



dC 
dT 



i\r'^^',c]~(w"^^^c+c 



■)+K"dge^ (B57) 
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Figure 14. Thermal leads: (a) Steady-state number of vibrons {n,^) 
along the ion chain in the absence of the static spin-vibron coupling 
AcOfj = (see the text for the remaining parameters). The bulk of the 
chain yields a homogeneous number of vibrons. The bars (yellow, 
red, and blue) correspond to the numerical solution of Eqs. (B57)- 
(B58). (b) Same as above but setting a strong spin-vibron coupling 
AcOfj = 200/; ^2„- This turns the homogeneous number of vibrons 
into a step-like function 



where the matrices W'^'^^'^ , K'^''^'^ have been defined in 
Eqs. (B 13). Due to the on-site energy shifts, we have to mod- 
ify the J matrices as follows 



dredge _ ~ cj 



J, 



'ajp ^ 



a,l5e{a,K,T}, (B58) 



In Fig. 14, we represent the mean value of vibrons in 
the steady state of a A^ = 13 ion chain, where we re- 
call that the chosen species are T =^"*Mg+, a =^^Mg+, 
and K =^Be+. We consider the following trap frequencies 
{(0ax,(0ay,(0az)/^'^ = (5,5,0.1)MHz, and the laser-cooling 
parameters Ai^ = — O.SFt, Ai3^ = —O.GTx, and i^L™ = 

It 

Hlsw — IAFt:, such that we expect each reservoir to thermal- 

13t 

ize to ni^ = 1.84, and ni3^ — 1.65. In Fig. 14(a), we represent 
our results in the absence of the on-site energy shifts Aco^ = 0. 
In analogy to the TQW [Fig. 1 1(a)], we recover the expected 
homogeneous mean number of vibrons along the whole bulk. 
In Fig. 14(b), we study the consequences of switching a very 
strong spin-vibron coupling A©,^ = 2007ij2cr- I" this case, 
the left half of the chain thermalizes to the left reservoir 



to the TQD described by ^lkr in Eq. (B53), the tunnel- 
ing of vibrons across the TQD also becomes rapidly rotat- 
ing in the regime of strong couplings Aco^ = 2007i^2<,' such 
that the current through the TQD is inhibited. In fact, we 
find numerically that the vibron current through the JC-ion is 
(/"''_^)ss ~ 1.6- lO^'^vibrons/i. This must be contrasted to 
the case of A(0^ = 0, where {Ip'^^)s,s, ~ 15vibrons/i. For ex- 
perimental time-scales, we can consider that the strong driv- 
ings Aco^ suppress completely the vibron current through the 
TQD. Hence, only the long-range tunnelings to the reservoirs 
influence the thermalization of the TQD ^lk-r ~ ^(c- 

//) Single-spin heat switch- We now describe a mechanism 
to switch on the vibron current across the TQD. We make use 
of the last ingredient in our toolbox (4), a periodic spin-vibron 
coupling (3) applied to the fc-ions 

H^^{t) ^ i(Aco+ + Ao)^ cJ,^Jcos(v^f - ^/c)";.,- (B59) 

According to Sec. A 4, and the explicit relations in Eqs. (A47)- 
(A48), we can achieve such a spin-vibron coupling by using 
a pair of laser beams with different frequencies (0\^ ^ ©2^, 
such that the frequency of the spin-vibron coupling (3) v^- = 
COliw = Oi^ — 0)2^. Moreover, by adjusting the laser intensi- 
ties, detunings, polarizations, and phases, we impose 



^A«„, 



^K = 0, Ao)^ = rA(0^ 



(B60) 



The idea is to use this periodic modulation to bridge the gra- 
dient of on-site energies between the two halves of the chain, 
assisting in this way the tunneling through the TQD. In the 
spin-independent case, the tunneling of vibrons can mimic a 
synthetic external gauge field [18]. In the current context, we 
exploit the spin-dependent drivings, such that depending on 
the parameter r, we can build a single-spin heat switch. 

Let us supplement the Liouvillian in Eq. (B53) with the pe- 
riodic spin-vibron coupling 



HLKK{t)^HLKK{t)=Hi^KR{t)+K{t)- 



(B61) 



In order to understand its effects, we move into another in- 
teraction picture with respect to the periodic driving Op^ -> 
Usyap^U^^, where U^v = e\p{ifQdt'H^^{t')}. This leads to 



,-if^(l+r(T;; )sin(v^f) 



, c. 



A< 

2v„ ' 



(B62) 



which can be inserted in the the tunneling of vibrons between 
the TQD and the leads (B54). By using the Jacobi-Anger ex- 
pansion for the first-kind Bessel functions ^n{z), namely 



^U sin 6 



LUzW" 



(B63) 



neZ 



together with the constraints (B60), it is possible to derive an 
effective Hamiltonian for the coupling of the leads to the TQD 



\"I(t/SS 

to («,„ 



hl for ia < 7k:, whereas the right half thermalizes -f^LicR *~ H '^JtaPK-SiyCKi^ '^''^pJ)'-'ia'^PK 



~ wr for ia > 7k-. We can thus conclude that the 
prediction for the behavior of the leads (B56) describes con- 
siderably well the actual steady state of the mixed ion chain. 
Let us also note that, according to the coupling of the leads 



'(!<Pk 



+ L 24p.ai(C^(l+m;J)flj'^flp^+H.c. 

ia>PK 



(B64) 
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where we have considered that all species have the same trap 
frequencies, and we have also used a rotating wave approx- 
imation for Ji^p^ ^ jjAOj^l. As announced previously, the 
expression above (B64) shows that for the resonance condi- 
tion Vk: = jAa)(7, the periodic spin-vibron coupling is capable 
of assisting the tunneling of vibrons across the TQD. More- 
over, the spin-dependence of the effective tunneling strengths 
via 5i (Ck-(1 + >''^pk)) ^^^ ^^ exploited to build the aforemen- 
tioned single-spin heat switch. By setting r = 1, we obtain 
3i(C^(l + (T^J) ==ai(2C(c)|t(c)(t)c|, which means that tun- 
neling is only allowed if the fc-ion is in the spin-up state. 
Therefore, by controlling the (C-spin using microwave or laser 
radiation (i.e. n pulses), it is possible to switch on/off the heat 
current. 

In order to check these predictions numerically, we consider 
a simplified setup, namely a a — K— a junction mimicking 
the connection of the thermal leads to the TQD. Rather than 
studying the steady state, we will concentrate on the coherent 
dynamics to show that the tunneling can be switched on/off 
by controlling the spin state of the the ic-ion. Let us define the 
parameters for this setup. We consider the usual trap frequen- 
cies (a)Q;x:fi^c(vift'c(,-)/27r = (5,5,0.1)MHz, which determine 
the equilibrium positions, and thus the tight-binding Hamilto- 
nian for the vibrons //tb. The static spin-vibron coupling is 






(B65) 



where Aa)+ = -2{J2^2^-J\„\^) and Ao)^ = 10^7i„2^- This 
provides an energy gradient that inhibits the tunneling across 
the TQD. The periodic spin-vibron coupling of the K"-ion 
//sy(f) is given by Eqs. (B59) and (B60), where we set r = 1. 
All these ingredients contribute to the dynamics given by 
H{t) = H" +//j!^(f) +//tb, which is solved numerically and 
compared to the theoretical predictions from H^^ (B64). 

We consider the initial state p(0) = Pi^ <8)p2^(8)P3(,, where 
Pia = l";a)("!al ® \^ia){^'a\ i^ determined by the vibrational 
Fock states rii^ = 1, n2,c ~ 0, nj^ — 0, and the spin states 
Sj^ G {tia'iia}- ^^ Want to Understand how the dynamics of 
such an initial state for si^ =iicr> ^3a =ti(T> is modified by 
changing S2^ G {12^,4-2^}- In Fig. 15(a), we set ^2^ =t2^, and 
observe how the vibron initially at the leftmost a-ion tunnels 
through the TQD until it reaches the rightmost d-ion. The 
solid lines correspond to the numerical solution of the full 
Hamiltonian H{t), while the open symbols are the predictions 
given by HI^. The agreement between both results supports 
the validity of our derivation. Therefore, we expect that by 
interspersing n pulses that invert the ic-spin |t2 ) ^ Ii2 )> we 
can switch on/off the vibron current. In Fig. 15(b), we show 
that two consecutive n pulses allow us to switch off the vibron 
current momentarily, which thus confirms our prediction. 



Appendix C: Spin-based measurements for heat transport 

The objective of this section is to present a detailed deriva- 
tion, supported by numerical simulations, of the Ramsey 
probes for measuring arbitrary vibronic observables (8). Then 
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Figure 15. Single-spin lieat switcli: (a) Mean number of vibrons 
(n,-^ ) as a function of time in the regime of photon-assisted tunnel- 
ing i'2^. =t2^ (see the text for the remaining parameters). The solid 
lines ((/ii„) blue, (n2^-) rd, («3^) yellow) represent the exact solu- 
tion of H{t), while the open symbols ((ni„) circles, («2r) diamonds, 
(n3^) squares) correspond to the effective photon-assisted-tunneling 
Hamiltonian ^lkR- ("'^ Mean number of vibrons («/^) (same as in 
(a)) as a function of time, where the KT-spin undergoes two consec- 
utive ;r-pulses that switch off/on the current. Note that the ;r-pulses 
are synchronized with the period of the spin-vibron coupling H^y{t). 



we particularize to the measurements of the vibron number 
and the heat current. 



1. Ramsey-type measurement of vibronic observables 

Let us start from the bulk spin-vibron model in Eq. (5), and 
consider a generic spin-vibron coupling for the fc-spins 



i/r(0^/?s''v=I^AoO,,CT4, 



(CI) 



where the "tildes" refer to the interaction picture with respect 
to the spin and on-site vibron Hamiltonians Hq = H" + H^ + 
//rvo- Here, we have introduced an arbitrary vibronic operator 
Oi^ ~ Oi^{{a- ,a- }), and the spin-vibron coupling ?^). In 
the sections bellow, we will specify to measurements of the 



vibron numbers (9;^ = n,-^. 



and vibron currents (9,v = /■ 



vib 



Since we want to probe the steady state of the bulk ion 
chain, the above spin-vibron coupling should disturb mini- 
mally the dynamics of the vibrons. Therefore, we impose 



l^ol«4 



■ A± 



(C2) 
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This allows us to partition the bulk dissipative spin-vibron 
model (5) into two terms 



-Sfo(Abulk) — — i[^rvt, Abulk] + ^bulk(Abulk), 
ibulk) 



-Sfl(Abulk) — — i[^sv'Abulk] 



(C3) 



where H^^i is tunneling part of the renormalized tight-binding 
model (5). The idea now is to project onto the steady-state 
of the bulk ions, which is given by -2o(/X^^[]^) = 0. We use 
again the projection-operator techniques in Eq. (A2 1 ), where 
the projector is now ^buik{»} = Mbuik^Tra.spinlTrvibd*}}. 
This allows us to obtain an effective master equation for the 
fc-spins, which reads as follows 



~ spin 



dA^' 



df 



■ °-^ Ramsey VMk" 



-spins. 



i[//R,An+i^R(Ar). (C4) 



Here, we have introduced a Hamiltonian that is responsible 
for the coherent part of the probe 



//R = ^iAo(0,Jss(T^ 



(C5) 



This term maps the information about the mean value of the 
vibronic operator <9;^ in the steady-state of the bulk ions onto 
the phase evolution of the spins. The vibronic fluctuations will 
be coded into the incoherent part of the probe 

^r(-)= L ilo%,,0,.(0)(c74.C7j^-C7^^C74.)+H.C.. 
IkJk 

(C6) 
Here, we have introduced the spectral function of the correla- 
tor between two vibronic observables 



So.,o,Sco)= df(O,,(f)O,,(0))sse 



+icot 



(C7) 



where the operators Oi^ = Oi^ — {Oi^)ss quantify the fluctua- 
tions from the steady-state values, and we use 
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Figure 16. Spin-based measurements for heat transport: (a) 

(upper panel) We consider a situation analogous to a thermal quan- 
tum wire (TQW) (i.e. a bar connected to two heat reservoirs with 
different temperatures), and try to find a probe of the temperature 
along the TQW. (mid panel) In the trapped-ion scheme [Fig. 1], we 
switch on the lasers L*,^ for the ic-species leading to a static and weak 
spin-vibron coupling (3). If the KT-spins are initialized in a linear su- 
perposition 1+,-^) = (Itiv) + \iiic))/V^, the spin dynamics resembles 
a Ramsey interferometer capable of capturing the information about 
the mean vibron number and its fluctuations (lower panel), (b) (up- 
per panel) We consider a situation analogous to a thermal quantum 
dot (TQD) connected to two thermal leads in equilibrium with two 
reservoirs held at different temperatures. We want to find a probe of 
the heat current flowing through the TQD. (mid pannel) In addition 
to the static spin-vibron coupling of the c-ions of (a), the lasers L'^ 
should induce now a periodic and weak spin-vibron coupling (3). In 
this case the driving is responsible for assisting the tunneling, but 
also for mapping the information about the vibron current to the spin 
coherences in a Ramsey-type interferometer (lower pannel). 



The second pulse, and the projective measurement, are equiv- 



((5;^(f)(5;^(0))ss=Tr{(9,„e-^o'(9;,,jUbuik}- (C8) alent to the measurement of the spin coherences {afXt)) 



Therefore, the zero-frequency component of the above spec- 
tral function (C7) determines the dephasing dynamics (C6) 
of the spins forming our probe due to the vibron fluctuations 
around the steady-state. Thus, by measuring the spin dephas- 
ing, we will acquire information about the quantum noise of 
the vibronic observable being probed. 

We now describe in detail how the mean value and the fluc- 
tuations of the vibronic operator Oj^ can be measured in anal- 
ogy to a Ramsey interferometer [Fig. 16(a)-(b)]. Let us an- 
alyze the case where the probe is made of a single JC-ion, 
which is initialized by a 7r/2-pulse in the spin state /^^^^'"(O) — 

H-,-J(+,J, where |+,J = i\%) + |i,,))/V2. Then, the bulk 
ions evolve under the Liouvillian J^q (C3), such that their 
vibrons reach the steady state, while the JC-spins evolve ac- 
cording to ^Ramsey (C4), acquiring thus information about 
the vibron observable 0,^. In order to recover this informa- 
tion, we perform another ;r/2-pulse, and measure the prob- 
ability of observing the K"-ion in the spin-down state P|. . 



where af^ = |t,-,> (i,,- |e-"«o ' + |;,. J (t, Je+'«'o '. According to 
Eq. (C4), the ic-spin coherences evolve as follows 



(af^(f)>=cos(Ao(0,JssOe 



-A2Re{So,^o,J0)}r 



(C9) 



Therefore, by measuring the spin populations as a function 
of time, we expect to get damped oscillations, the period 
of which gives us information about the mean number of 
vibrons, while their damping is proportional to the vibron- 
number fluctuations in the steady-state. Let us note that 
the spin-population measurements can be performed through 
the state-dependent fluorescence of the trapped ion, a tech- 
nique routinely used in many laboratories that allow for ac- 
curacies reaching 100% for detection times in the millisec- 
ond range [19]. Let us remark that, since we are interested 
in steady-state properties of the vibrons, this measurement 
scheme is not sensitive to the time-resolution of the spin-state 
readout. Hence, this does not pose any limitation to the target 
accuracies reaching 100%. 
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. Let us finally note that, according to Eq. {CI), if the probe 
consists of several K"-ions, we will also have access to the two- 
point correlations of distant ions. 



2. Particular applications: vibron number and current 

a. Measurement of the vibron number 

In order to tailor the coupling (CI ) to probe the vibron num- 
ber (i.e. Oi^ ~ tii^), we must resort to a weak and static spin- 
vibron coupling (3). According to Eq. (B31), we can achieve 
such a spin-vibron coupling by using a pair of laser beams 
with equal frequencies, leading to 



K ^Y.\^^Kn,^(^l 



(CIO) 



where we recall that the Rabi frequency and the Lamb- 
Dicke parameter, defined by Eqs. (A42)-(A43), can be fully 
controlled experimentally. In light of the notation used in 
Eq. (CI), we identify Oj^ — tii^ [Fig. 16(a)], and Aq — Ao)^, 
which can be tuned to fulfill the required probe condition (C2). 
If we restrict to a single probing ion labeled by p^, according 
to Eqs. (C4), the coherences evolve as follows 



(^P.(0)=cos(Aa)^(npJssOe 



-(Ao^)2Re{S„p^„p^(0)}; 



(Cll) 



which coincides with the description in the main text, and al- 
lows us to extract the quasi-particle mean number and fluc- 
tuations. We remark that for more than one probing ion, the 
dissipative dynamics would also gives us information about 
the vibronic correlations between distant ions. 

In order to support our derivations, let us analyze numeri- 
cally the Ramsey measurement for the TQW. Due to the intro- 
duction of the fc-spins, the dynamics of the system is no longer 
quadratic as in Sec. (B 2), which forbids finding a closed sys- 
tem of differential equations for the vibronic two-point cor- 
relators. Therefore, we have to obtain numerically the time 
evolution of the complete density matrix /Xbuik(0 given by 
Eq. (C3), and then calculate the observable {of (f)). Due to 
the computational cost of this problem, let us simplify max- 
imally the setup where the Ramsey measurement can be de- 
veloped by considering a thermal quantum dot (TQD) once 
more. However, in contrast to Sec. B 2 a, we will consider 
the arrangement x —K— x, where x =^^Mg+ and K =^Be+. 
We consider the same parameters introduced in previous sec- 
tions, but set the detunings Aj^ = — 0.6rT, At,^ = — O.Sr^, and 
the Rabi frequencies i^Lsw = ii^sw = F^ for the laser cooling 

It 3t 

of the T-ions, where we recall that Tt:/2k = 41.4MHz. With 
these parameters, the effective cooling rates (A28) of the X- 
ions would be 7ij/2;r?a 86kHz, a.nd'Yj,j2n^ 106kHz. Addi- 
tionally, the mean number of vibrons for each reservoir (A29) 
would be h\^ = 1.65, and ns^ = 1.63. The trap frequencies 
are (a)a.r,«av,fi^te)/2?r = (5,5,0.25)MHz, which lead to a 
vibron tunneling strength of Ji^^^l'^'^ — -^3x2^/2?^ ~ 35kHz. 
With these parameters, we expect that the T-ions thermalize 
very fast and act as an effective reservoir for the bulk )C-ion. 
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Figure 17. Ramsey-type measurement of the vibron number: (a) 

Dynamics for the coherence of the probe spin {a\ (t)) in a Ramsey- 
type measurement. The grey dashed line would represent the pe- 
riodic oscillations in Eq. (C22) that would be obtained in a noise- 
less scenario. However, due to the quantum noise S„; „, (0), such 
oscillations get damped as shown by the numerical solution (green 
solid line). The crosses correspond to a numerical fit {a^ if]) = 
cos(a?)exp{— to}, with fitting parameters a.b, which allow us to re- 
cover the mean value and fluctuations of the number of vibrons in 
the steady state via Eq. (C22). (b) Mean value of the vibron num- 
ber obtained from the numerical fit ("2k)ss = fl/Aft)^ (solid line) 
as a function of the probing strength. The dashed line represents 
the theoretical prediction (CI 4). As expected, for Aa)^^ — > 0, the 
probe does not disturb the properties of the bulk vibrons, and we re- 
cover the predicted mean number of vibrons (CI 4) (dashed line), (c) 
Quantum noise of the vibron number obtained from the numerical fit 
(S,!2 ni (0))ss = ^/(Aco^) (solid line) as a function of the probing 
strength. For Aco^^ — ?> 0, the probe does not disturb the properties of 
the bulk vibrons, and we recover the prediction (C14) (dashed line). 



According to Eq. (5), the K"-ion evolves under /ibuik 
-S?2<,2<,(j"buik), where 



■^2.2, (.) 






K' ^K ^K ' 



(C12) 
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where the dissipator contains the heating and coohng terms 



bulk 



(•) =A^ 



K^K \ -^K '-■K '^K '^K > 
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Let us note that, in this case, the mean number of vibrons 
and the noise fluctuations in the steady state can be calculated 
exactly by means of the quantum regression theorem 
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(C14) 
Therefore, this particular TQD offers a neat playground to test 
the proposed measurement scheme. 

We now solve numerically the master equation (CI 2) con- 
sidering the above realistic parameters, and compute the dy- 
namics of the coherences (ffj (f)). In Fig. 17(a), we repre- 
sent the numerical results for these coherences (solid line), 
and perform a numerical fit (crosses) to the expected behav- 
ior in Eq. (C22), which allows us to infer the mean value and 
the quantum noise of the vibron number {n2i^)ss, 5'„, „2 (0). 
In Figs. 17(b)-(c), we represent the results obtained from 
this numerical fit as a function of the probing strength Aco^ . 
When the probe is too strong, there is an important back- 
action on the bulk vibrons, and the values obtained from the 
fit depart from the theoretical prediction (CI 4). Conversely, 
for Aco^ -^ 0, the probe disturbs minimally the system, and 
one attains an acceptable agreement with the theoretical pre- 
dictions (CI 4). It is important to emphasize that, although 
Aco^ — 7> 0, the time required for the Ramsey measurement is 
fixed to f = 10ms, which is a reasonable regime considering 
typical decoherence times in trapped-ion experiments. Hence, 
the limit A©^ — > does not require to maintain the spin co- 
herences for prohibitively large experimental times. 



h. Measurement of the vibron current 

We now address a possible way of designing the cou- 
pling (CI) to probe the vibron current (i.e. (9,-^ = ^T\^- ^^ 
particular, we analyze the current through a TQD connected 
to the reservoirs by a couple of leads [Fig. 16(b)]. Let us re- 
call that this setup is described by the Liouvillian in Eq. (B49), 
where the reservoirs correspond to the laser-cooled T-ions, the 
leads to the CJ-chains, and the TQD to the fc-ion. As dis- 
cussed below Eq. (B49), the current through the TQD can be 
suppressed by means of a strong energy off-set between the 
two halves of the chain. Then, a periodic spin-vibron cou- 
pling (B59) serves as gadget to switch on the current, such 
that the tunneling strengths depend on the particular spin state 
of the K-ion (see Eq. (B64)). 

In this section, we will make use of this spin-dependence 
to build a Ramsey probe for the vibron current. However, in 
order to design design the Ramsey-type coupling (CI) such 
that Oi^ = I]*Lf, we need to modify the spin-vibron couplings. 



In particular, we exploit a bi-chromatic spin-vibron coupling 
^sv(0= L 5(A<«+A«^.«f^/'Jcos(v^^„f-9^,„)np,. 

n=1.2 

(C15) 
Moreover, by adjusting the laser intensities, detunings, polar- 
izations, and phases, we impose 
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(C16) 



where we recall that Aco^, is the energy off-set between the 
two halves of the chain (i.e. the leads). 

In analogy to the derivation of Eq. (B64), to understand 
the effects of the bi-chromatic spin-vibron coupling (CI 5), we 
move into an interaction picture with respect to the driving 

a,,^ -^ «P. = -flp^e-'f-i^'"(^'^."-''/2)g-i?,.2a,5^sm(v^.20, 

(C17) 
By using the Jacobi-Anger expansion again, together with the 
constraints (CI 6), it is possible to derive an effective Hamil- 
tonian for the coupling of the leads to the TQD 
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(CI8) 



where we have used a rotating wave approximation for 
JiaPK ^ jl^^cT I- As announced previously, the expression 
above shows that for the resonance conditions Vi^i = V^^ ~ 
^AcOfj, the bi-chromatic spin-vibron coupling is capable of 
assisting the tunneling of vibrons across the TQD. 

We will now exploit the spin-dependence of the effective 
tunneling strengths via the Bessel function 3m(CK:,20'L) to 
build a Ramsey probe of the vibron current. In particular, 
taking into account that ^^2 <C 1, we can rewrite //lk:r — 
^L^R +^sv' where we have introduced the Hamiltonian 

(C19) 
This term describes a spin-independent tunneling of vibrons 
across the TQD, which will be responsible for setting a vi- 
bron current. The important feature of the assisted-tunneling 
strength is that it has become purely imaginary, which be- 
comes relevant in the definition of the current operators (B21). 
This turns out to be crucial to devise the Ramsey probe, since 
the remaining terms in the Hamiltonian can be written as 

^ = 3A/(/;^+rj^Ja^^, (C20) 

where we have introduced the dimensionless coupUng 

2(5o(7r)+a2(7r)) 
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Figure 18. Effective tunneling for tlie current Ramsey measure- 
ments: (a) Mean number of vibrons (n,-^) as a function of time in 
the regime of photon-assisted tunnehng S2^ =^1^ (see the text for the 
remaining parameters). The solid lines ((ni„) blue, («2i-) red, («3„) 
yellow) represent the exact solution of H{t), while the open symbols 
((;ii^) circles, («2i-) diamonds, (ny,^) squares) correspond to the ef- 
fective photon-assisted-tunneling Hamiltonian //lktR = ^lk-r + ^sv 
(b) Time evolution of the spin coherence of the probe K-ion, (cTj (t)}, 
for an initial spin state ll'aKC^)) = 1+2^)- Due to the effective spin- 
current coupling (C20), the the coherences display periodic oscil- 
lations that depend on the periodicity of the assisted tunneling of 
vibrons between the ions. 



Remarkably enough, we can make the coupling of the Ram- 
sey probe arbitrarily small by simply letting ^^.2 -^ 0. In this 
limit, we have an ideal Ramsey probe of the vibronic current. 
According to Eqs. (C4), the fc-spin coherences evolve as 



(a^^(O) =cos(l/(n';,^),s,Oe-^''''^^''''^'''^^''»', (C22) 



where we have made use of the fact that (/"'j!,^)ss = {^p'^^)ss 
in the steady state. We have also defined the zero-frequency 
intensity noise 



S'P.'J^) -= l^dt{IpAt)tpAO)U, (C23) 



where we have introduced the current fluctuations through 
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Let us note that, in order to give supporting numerical evi- 
dence of this prediction, the minimal setup to explore would 
be a T — (7 — ic— (7 — T chain. In analogy to the single-spin 
switch, computing the dynamics of the spin coherences in this 
case becomes a non-linear problem that exceeds our numer- 
ical capabilities. Therefore, we cannot obtain the analogue 
of Fig. 17 for the current operator. Instead, we will con- 
tent ourselves with showing that the effective Hamiltonians in 
Eqs. (C19)-(C20) describe the dynamics of a C7 — K" — <T setup 
faithfully. This has been shown in Fig. 18. 
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